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PREFACE 



The PDP-8 comes to the user complete with an extensive selection of system programs 
and routines making the full data processing capability of the new computer immediately 
available to each user, eliminating many commonly experienced initial programming 
delays. 

The programs described in these abstracts come from two sources, past programming 
effort on the PDP-5 computer, and present and continuing programming effort on the 
PDP-8. Thus the PDP-8 programming system takes advantage of the many man-years 
of program development and field testing by PDP-5 users. 

Although in many cases PDP-8 programs originated as PDP-5 programs, all utility and 
functional program documentation is issued in a new, recursive format introduced with 
the PDP-8. 

Programs written by users of either the PDP-5 or the PDP-8 and submitted to the users' 
library (DECUS - Digital Equipment Corporation Users' Society) are immediately avail- 
able to PDP-8 users. 

Consequently, users of either computer can take immediate advantage of the continuing 
program developments for the other. 
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CHAPTER 1 



INTRODUCTION 



The PDP-8 Floating-Point System enables the programmer to concentrate on the logic of his computation 
rather than on decimal points. The system maintains a constant number of significant digits throughout 
the computation, thereby enhancing the accuracy of the result. 

Floating-point notation is particularly useful for computations involving numerous multiplications and 
divisions where magnitudes are likely to vary widely and where only crude predictions can be made as to 
the amount of variation involved. The main advantage of the PDP-8 Floating-Point System is derived 
from the ability to store very large or very small numbers by storing only the significant digits together 
with the exponent for that number. In an integral or fractional fixed-point machine, the programmer 
must either include all of the O's between the significant digits and the decimal point (and in so doing 
loose considerable accuracy), or account for the point in each step of his programming. 

The PDP-8 Floating-Point System is constructed as a self-contained package which includes its own input, 
arithmetic, and output routines. It allows the programmer to use floating-point arithmetic without having 
to construct his own arithmetic subroutines. For example, the polynomial: 

Y = AX 3 + BX 2 + CX + D 

may be programmed as follows (the variables X and Y, and the parameters A, B, C, and D are stored in 
registers of the same name): 



JMS I 7 
FGET A 
FMPY X 
FADD B 
FMPY X 
FADD C 
FMPY X 
FADD D 
FPUT Y 
FEXT 



/ENTER FLOATING PACKAGE 

/LOAD PSEUDO AC 

/MULTIPLY 

/ADD 

/MULTIPLY 

/ADD 

/MULTIPLY 

/ADD 

/STORE 

/EXIT FLOATING PACKAGE 
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CHAPTER 2 
FLOATING-POINT REPRESENTATION 



A floating-point number consists of a mantissa and an exponent. For example, the decimal number 12 may 
be represented as: 

1200. 0-10" 2 

120.00-10" 1 

12.000-10° 

1.2000-10 1 

.12000-10 2 

etc. 
where the exponents are the numbers —2, —1, 0, etc. 

Since the PDP-8 is a binary machine, floating-point numbers are stored as floating-point binary internally. 
For example, the binary number 1 10 (6 decimal) may be represented as: 

11000.0-2" 2 (24-1/4 = 6) 

1100.00-2" 1 (12-1/2 = 6) 

110.000-2° (6-1 =6) 

11.0000-2 1 (3-2 = 6) 

1.10000-2 2 (3/2-4 = 6) 

.110000-2 3 (3/4-8 = 6) 

etc. 

Notice that the binary exponent is always a signed integer. The PDP-8 Floating-Point System uses the 
3> convention:^ fe ^(^%^ I 

^- \j 1/2 < |mantissa| <i. \\'-Jl tM ^CLi*»\£ib.M+>i {~~) °*fy) t8^4L~*4***AZ^fL*4~* 

1\ When this is true, the word is said to be normalized. *■ l ^/ 

The value of the number is then: 

MANT.SSA-2 EXPONENT 

where the MANTISSA is a signed quantity. The result of this convention is that more significant bits are 
retained. For example, the number .10 (decimal) is equal to: 
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0.00 Oil 001 100 I 110 Oil 001 100 I no on obi 100 110 



in bin 



ary . 



O- 0&S}4CS\4-6, 3/4- 



When this word is stored in two 12-bit words (24 bits), the leading 0's a»e nonsignificant and only 20 bits 
of significance are maintained. If the number is rewritten as: 



0.11 001 100 110 



011 001 100 110 



011 001 100 110 .... 



and the mantissa stored in two 12-bit words, 23 bits of significance are maintained. 
The exponent is stored in a third word making a total of three words for storage. 



sign of exponent 













sign o 


" mantissa 





















binary exponent 2's complement 
signed quantity 



high order mantissa word 
low order mantissa word 
The number .1 (decimal) would be stored as: 

111 111 111 101 | 011001100110 | 011001100110 
or written in octal: 

7775 3146 3146 

The address of this number would be considered 2146 if it were stored as follows: 

2146/7775 
2147/3146 
2150/3146 

6 3 / a- C 3 ' 4- 
The number -.1 (decimal) would be -.063 1463 14 (octal) which is 0.000 110 011 001 100 110 011 001 100 

(binary). 



-V0l„ Normalized, this becomes: (2~°) '(0.11 001 100 110 011 001 100 1 10) 2 
Complemented and incremented to become the negative 2's complement 



M * 'o this becomes: 

ro 

which is: 



(2 ) -(1.00 110 011 001 100 110 011 010) 2 
(2" 3 ) 8 -(4 6 3 1 4 6 3 2 ) 8 
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This number would be stored (within octal) as: 
7775 exponent of -3 

4631 high order mantissa word 

4632 low order mantisaword 

Since the mantissa is greater (in magnitude) than 1/2 and less than 1, its binary point is considered to 
be between bit and bit 1 of the high order mantissa word. 

Further examples: 

12 (decimal) 

1100. = (2 4 ) (0.11,000,000000000000000000) 

OOOCi- y O 'J o 

would be stored as (written in octal): 



0004 










3000 










0000 










Octal 


Decimal 




Floating Binary 




-1 
-10 


-1 
-8 


S 0001 

C OOOo 

0004 


6000 
6000 


0000 
0000 


.4 


.5 


0000 


2000 


0000 


3 


3 


0002 


3000 


0000 


14 


12 


0004 


3000 


0000 


.2 


.25 


7777 


2000 


0000 


14.26074475 


12.3456 


0004 


3054 


171Z 


-14.26074475 


12.3456 


0004 


4723 


606$ 



? fJo. 



Convert the decimal number 1 .6 to PDP-8 floating-binary form for storage. 
An easy method to convert a decimal fraction to an octal fraction is illustrated below: 

1.6 1Q = 1.46314631 

= i.ioGJiiojoii|oo1ioqiioioiijooi 
= (2 1 )(o.i^ooi|ioandoi looi.ioQi io) 

= 0001 
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X 

4) 


.6 

8 

.8 


X 

6) 


8 

.4 


X 

3) 


8 
.2 





X 


8 


1) 


.6 




8 


4) 


.8 




8 


6) 


.4 




8 


3) 


.2 




8 


1) 


.6 



3146 
3146 



ARITHMETIC 

Since floating-point numbers are stored in a three-register format, the floating-point system uses a "psuedo" 
floating accumulator (FAC) which consists of three registers in the floating-point package: 44, 45, and 46. 
Register 44 contains the exponent; 45 and 46 contain the high and low order parts of the mantissa, respec- 
tively. 

BASIC FLOAT I NG-POINT COMMANDS 

The basic floating-point commands include the following: 

load floating accumulator 
store floating accumulator 
add to floating accumulator 
subtract from floating accumulator 
multiply by floating accumulator 
divide into floating accumulator 
normalize floating accumulator 

All arithmetic operations are called through an interpreter. The command codes have a format that is 
almost identical to the format of the PDP-8 memory reference instruction, namely: 



OP 


1 


z 


ADDRESS 



where the op code is from 000 to 1 1 1 or o - 7 . 

o o 

INTERPRETER 

The interpreter contains, at all times, the address of the memory location containing the next pseudo in- 
struction to be executed. This is initially stored when the program enters the interpreter using a JMS 
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instruction. When the interpreter encounters an instruction with an op code of and with bits 8-1 1 of 
the pseudo instruction equal to 0, it exits to the next memory location. 

FLOATING-POINT INSTRUCTIONS 
The floating-point instructions are: 

Op Code Mnemonic Effect 

1 FADD Floating Addition 

Add the contents of the effective address 
to the floating accumulator. 

2 FSUB Floating Subtract 

Subtract the contents of the effect address 
from the floating accumulator. 

3 FMPY Floating Multiply 

Multiply the floating accumulator by the 
contents of the effective address. 

4 FDIV Floating Divide 

Divide the floating accumulator by the 
contents of the effecitve address. 

5 FGET Floating Get 

Load the floating accumulator with the 
contents of the effective address. 

6 FPUT Floating Put 

Store the contents of the floating accumulator 
at the locations specified by the effective 
address. The contents of the floating accum- 
ulator are unchanged. 

7 FNOR Floating Normalize 

Normalize the contents of the floating 
accumulator. 

is decoded as follows: 

Bits 8-1 1 = 0000 Floating Exit 

Return control to following instruction 

= 0001 Floating Square 

Square the contents of the floating ac- 
cumulator 
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= 0010 



= 0011-1111 



Floating Square Root 

Take the root of the absolute value of the 

floating accumulator 

Expandable commands 



This may be summarized: 

FADD Y; 1000, 

FSUB Y; 2000, 

FMPY Y; 3000, 

FDIV Y; 4000, 

FGET Y; 5000, 

FPUT Y; 6000, 
FN OR ; 7000, 
FEXT ; 0000, 
SQUARE ; 0001 
SQROOT; 0002 



C(FAC) + C(Y) -> C(FAC) A 

C(FAC) - C(Y) -> C(FAC) I 

C(FAC) x C(Y) -> C(FAC) f 

C(FAC) t C(Y) -> C(FAC) J 
C(Y) -> C(FAC) 



Result is 
normalized 



C(Y) 
unchanged 



C(FAC) -> C(Y) 

C(FAC) normalized -> C(FAC) 

exit from interpreter to instruction following this command 

C(FAC) 2 -> C(FAC) 

C(FAC) 1/2 _> C(FAC) 



The assembler recognizes all of these mnemonics except for SQUARE and SQROOT which may be defined 
as: 

SQUARE=0001 
SQROOT=0002 

The floating-point interpreter is entered with an effective JMS 5600. 



For example: 








SQROOT=0002 






*7 








5600 








*200 










JMS 1 7 


200/ 


4407 




FGET A 


201/ 


5206 




SQROOT 


202/ 


0002 




FPUT 1 B 


203/ 


6611 




FEXT 


204/ 


0000 




HLT 


205/ 


7402 


A, 


0003 


206/ 


0003 




2000 


207/ 


2000 




0000 


210/ 


0000 


B, 


300 


211/ 


0300 


$ 
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When this program is started at 0200, it will halt at location 205. The state of the machine will be: 

44/ 0002 

45/ 2000 floating accumulator contains 

46/ 0000 2.0, i.e., (4.0) 1 /2 

206/ 0003 

207/ 2000 register A contains 4.0 



210/ 0000 

300/ 0002 

301/ 2000 

302/ 0000 



300/ 0002 

301/ 2000 answer stored here 
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CHAPTER 3 
FLOATING-POINT INPUT/OUTPUT 

FLOATING-POINT INPUT 

The basic floating-point package contains an input routine to read characters from the 33 ASR keyboard. 
Numbers are read in in decimal and converted to the internal floating-point binary format. Input format 
is floating decimal . The number 726.7 may be read in in any of the following forms: 

726.7 

.7267E3 

.7267E+03 

+7267E-1 

etc. 

Input is terminated when a character is typed that is not a part of this format. The conversion of "12.0." 
would be terminated on the second "." and the binary number: 

000000000100 (octal 0004) 

011000000000 (octal 3000) 

000000000000 (octal 0000) 

would be in the floating-point accumulator upon completion of the conversion. 

Flags 
There are several flags associated with the input routine that are useful to the programmer. 

0056 This register is a switch that has the following meaning: 

If C(56) = 0, do not type a line-feed after a carriage-return was read. 
If C(56) = 7777 , type a line-feed when a carriage-return is inputted. 
This switch is initially set to 7777 . 

0057 This register contains the character that terminated the input con- 
version. 

0060 This register contains 0000 if no input conversion was made; i.e., a 

space or other terminator was initially typed. 
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The input routine is entered with an effective JMS 7400. It returns control to the instruction following 
the calling JMS upon receipt of a terminator. The floating accumulator contains the input number in 
normalized floating-binary. Register 0057 contains the terminating character in ASCII, and C(0060) in- 
dicates whether or not there was a valid input. 

For example: 

*5 

7400 
*7 

5600 
*200 

JMS I 5 /INPUT ROUTINE 

JMS I 7 /CALL FLOATING POINT 

FPUT A 

FEXT 

JMS I 5 

JMS I 7 

FPUT B 

FEXT 

HLT 

A, 



B, 



$ 

When this program is started at 0200 and the following is typed: 
X2.0Y 

The program will halt at location 0210, and A and B will contain 

A, 


and C(57) = 0331 , the second TERMINATOR 

B, 0002 
2000 
0000 

The first input was considered a because a terminator "X" was inputted. 
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This program could be written to ignore the non-numeric information as follows: 



*5 




7400 




*7 




5600 




*200 




JMS 15 


/CALL INPUT ROUTINE 


TAD 60 


/ANY VALID INPUT? 


SNA CLA 




JMP .-3 


/NO - IGNORE 


JMS 1 7 


/YES 


FPUT A 


/STORE IT 


FEXT 




JMS 1 5 


/GET NEXT 


TAD 60 




SNA CLA 


/VALID? 


JMP .-3 


/NO - IGNORE 


JMS 1 7 


/YES 


FPUT B 


/STORE IT 


FEXT 




HLT /HALT 




A, 














B, 















Register 57 may be used for integrating control characters into the input. 

Rubout 

There is one special input character, the rubout. If it is struck before an input delimiter, the input rou- 
tine is restarted and previous numbers deleted. 

For example, assume the input routine has been entered. 

4- 

3 

276 -§1T 

ai. 

The input routine would exit with 1 (decimal) in the floating accumulator or: 

44/ 0001 

2000 
0000 
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FLOATING-POINT OUTPUT 

The basic floating-point package contains an output routine. Upon entry the contents of the floating ac- 
cumulator are converted to floating-point decimal and typed out on the ASR33in the following format: 

±O.XXXXXXXE±XX 

For example, if the floating accumulator contained: 

44/ 0002 

2000 
0000 

and the output routine were entered, 

-K).2000000E-K)1 

would be typed. 

Entry 

The output routine is entered by an effective JMS 7200. After outputting the contents of the floating ac- 
cumulator, it returns control to the instruction following the calling JMS instruction. The contents of the 
floating accumulator are lost. The floating-output routine has a flag or program switch in location 0055 
on page 0. 

If C(0055) does not equal 0, type a carriage-return/line-feed following each output. This location initi- 
ally contains 7777 . 

Program Example 

/SOLVE THE QUADRATIC EQUATION 

/AX 2 +BX+C=0 FOR VALUES 

/OF A, B, C. TYPE OUT ROOTS 

SQROOT=0002 

SQUARE=0001 

*5 7400 /POINTER TO INPUT ROUTINE 

7200 /POI NTER TO OUTPUT ROUTI NE 

5600 /POINTER TO INTERPRETER 

*200 

KCC 
TLS 

BEGIN, JMS I 5 /INPUT A 

JMS I 7 /ENTER INTERPRETER 
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FMPY TWO 


/2.0 




FPUT A 


/STORE 2.0- A 




FEXT 






JMS 15 


/INPUT B 




JMS 1 7 


/ENTER INTERPRETER 




FPUT B 


/STORE B 




FEXT 






JMS 15 


/INPUT C 




JMS 1 7 


/ENTER INTERPRETER 




FMPY TWO 


/MULTIPLY BY 2. 




FMPY A 


/MULTIPLY BY2A 




FPUT TEMP 


/4AC 




FGET B 






SQUARE 


/B SQUARED 




FSUB TEMP 


/B SQUARED -4AC 




SQROOT 


/TAKE SQUARE ROOT 




FPUT TEMP 


/STORE IT 




FGET B 






FMPY MINI 


/-1-B = -B 




FPUT B 


/STORE 




FSUB TEMP 


/-B-SQUARE ROOT 




FDIV A 






FEXT 


/ANSWER IN FAC 




JMS 1 6 


/OUTPUT IT 




JMS 1 7 


/ENTER INTERPRETER 




FGET B 


/-B 




FADD TEMP 


/-B + SQUARE ROOT 




FDIV A 






FEXT 


/ANSWER IN FAC 




JMS 1 6 


/OUTPUT IT 




JMP BEGIN 


/CONTINUE 


TWO, 


0002 
2000 
0000 


/CONSTANT 2.0 


MINI, 


0001 
6000 
0000 


/CONSTANT -1 .0 


TEMP, 









A, 









B, 









c, 








$ 








Both the input-conversion routine and the output-conversion routine use a typeout subroutine that contains 
the instructions: 
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TSF 

JMP .-1 
TLS 
CLA 

This means that the teleprinter flag must be set when these routines are entered or it must be in the process 
of typing. A TLS instruction can be put into the initialization section of a program. 

THE BASIC PACKAGE 



Subroutines 

When the floating-point interpreter encounters a op code, it further decodes bits 8-11 . If these bits 
equal 0, the interpreter exits. If the bits are nonzero, they are used in a table look-up to specify the 
address of a subroutine. The called subroutine may use the floating-point interpreter, the input, or the 
output, but it may not use bits 8-1 1 to call another subroutine. The look-up table is on page 16 of the 
interpreter listing. For example, if the user wished to define a routine that negated the contents of the 
floating accumulator, the following steps would be taken: 

On page 6 of the interpreter listing, there is a routine to negate the floating accumulator. It is a sub- 
routine and has, as entry point, the address 6000. If the negate subroutine were to be called by 0010 in 
the interpreter mode, the word 6000 would be placed in address 6554 of the calling table. Input and out- 
put could be called in the interpreter mode if 7400 were placed in 6555 (001 1) and 7200 were placed in 
6556(0012). 

NEGATE =0010 
INPUT =0011 
OUTPUT =0012 
SQUARE = 0001 
*7 

5600 
*200 

KCC 

TLS 

JMS I 7 /ENTER INTERPRETER 

I NPUT /CALL I NPUT ROUTI NE 

SQUARE /SQUARE IT 

NEGATE /CALL OUTPUT ROUTINE 

OUTPUT /CALL OUTPUT ROUTINE 

FEXT /EXIT 

HLT /HALT 
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DESCRIPTION OF BASIC FUNCTION 



Addition 



Floating-point addition is carried out by first aligning the binary points of the two numbers. This is ac- 
complished by scaling the smaller number to the right. Then the mantissas are both scaled right once so 
that overflow will not occur into the sign bit. A 2's complement addition of the mantissas is then made. 
The result is normalized and control returns to the interpreter. This may be represented as: 



Can operands be aligned? 



"1 

Yes No 



Align Exponents 

Add Mantissas Put larger number in FAC 

Normalize 

Return to Interpreter 4 



Subtraction 

Floating-point subtraction is accomplished by negating the operand, and then calling the addition sub- 
program. 

Negate Operand 
Call Addition 

Multiplication 

Floating-point multiplication is accomplished by adding the exponents together and then performing a 
double-precision multiplication. The result is normalized and control returns to the interpreter. 

Add Exponents 

Multiply Mantissas carrying result to 35 bits 

Normalize 

Return to Interpreter 

Division 

Floating-point division is accomplished by subtracting the exponent of the divisor from the exponent of 
the dividend. The mantissa is divided and the result is normalized. Control returns to the interpreter. 
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Is divisor = 0? 



sor = u f 1 

No Yes 

Subtract Exponents 

Divide Mantissas +1 to error flag 

Normalize Set result to highest possible 

Return to interpreter ^ 1 

Square 
The square routine calls the multiplication routine internally. 

Square Root 

The square root is calculated using Newton's method in which an initial approximation is made and then 
each succeeding approximation is calculated. The routine exits when two successive approximations are 
equal to within the least significant bit of the mantissas. 

Form first approximation, X , of v^NT 

then: 

X.+1 = 1/2 (X. + £L) 
i 



is I X. - X . | < least significant bit? 

i 

No 

Yes 



1 

r es 
Return to interpreter 

Error Flag 
Division by causes C(61) to be incremented by 1 . The quotient is set to the highest positive number. 

Attempting to extract the square of a negative number causes C(61) to be incremented by 1 . The root of 
the absolute value is taken. 

The contents of 61 are set to at the beginning of each divide and square root operation. 
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SUMMARY OF BASIC PACKAGE 



5600 
7400 
7200 



Entry Points 



Arithmetic Interpreter 
Floating Input 
Floating Output 



NOTE: Both the input and the output routines require that register 7 contains 
the interpreter entry point: 5600. 



55 
56 
57 
60 
61 



Flogs 

If ^ 0, type carriage-returry/line-feed after output. 

If ^ 0, follow each input carriage-return by a line-feed. 

Contains the input terminating character. 

Equals if no valid input. 

Is nonzero if (a) divide by or (b) square root of a negative number. 



FADD 

FSUB 

FMPY 

FDIV 

FGET 

FPUT 

FNOR 

FEXT 

SQUARE 

SQROOT 



1000 
2000 
3000 
4000 
5000 
6000 
7000 
0000 
0001 
0002 
0003 

0017 



Commands 

Floating Add 
Floating Subtract 
Floating Multiply 
Floating Divide 
Floating Get 
Floating Put 
Floating Normalize 
Floating Exit 
Square 
Square Root 

Expandable 
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Storage 
The floating accumulator is in registers: 

44 - Exponent 

45 - High Order Mantissa 

46 - Low Order Mantissa 

The basic floating-point package (Digital 8-5A-S) uses memory locations: 
7; 40-61; 5600-7577 

FIXED TO FLOATING/FLOATING TO FIXED 
Since the floating-point package stores numbers in the following format: 



Signed Exponent f 2's Complement Signed Mantisa 



T 2's Complement Sic 
Llmplied Binary Point 



and since the exponent indicates where the real exponent is, this information may be used to either con- 
vert a floating-point number to fixed point or to convert a fixed-point number to floating point. 

For example, assume that there is an integer in the accumulator that is less in magnitude than 2047. To 
float this number, the following sequence of steps may be employed: 



C, 



DCA45 


/PUT INTO HIGH-ORDER MANTISSA 


DCA46 


/PUT INTO LOW-ORDER MANTISSA 


TAD CI 3 


/ll (10) INTO 


DCA44 


/EXPONENT 


JMS 1 7 


/CALL INTERPRETER 


FNOR 


/NORMALIZE 


FEXT 


/LEAVE INTERPRETER 


0013 


/ll (DECIMAL) 



C13, 

At point C, we have set the binary point of the integer to the right end of the high order mantissa word, 
or eleven (decimal) locations to the right of the implicit binary point. 

To float this number, the floating accumulator is scaled right until the exponent contains the location of 
the desired binary point. To fix the floating accumulator as an 11-bit signed integer, the following se- 
quence of coding may be employed: 



3-10 



GO, 



DONE, 



CLA 

TAD 44 /FETCH EXPONENT 

SZA SMA /IS THE NUMBER <1 ? 

J MP .+3 /NO: 

CLA /YES: FIX IT TO 

JMP DONE+1 

TAD Ml 3 /NO: SET BINARY POINT AT 

SNA /l 1 (10) PLACES TO RIGHT OF CURRENT POINT 

JMP DONE /IT IS ALREADY THERE: ALL DONE 

SMA /TEST TO SEE IF IT IS TOO LARGE 

JMP ERROR /YES: NUMBER >2**11 

DCA44 /NO: SET SCALE COUNT 

CLL /0 TO C(L) 

TAD 45 /FETCH MANTISSA 

SPA /IS IT <0? 

CML /YES: PUT A I IN LEFT BIT 

RAR /SCALE RIGHT 

DCA45 /RESTORE IT 

ISZ 44 /TEST IF SHIFTED ENOUGH 

JMP 00 /NO: CONTINUE 

TAD 45 /ANSWER IN C(AC) 



Ml 3, -13 



/-ll (DECIMAL) 



This may be coded as a subroutine. 



EXTENDED FLOATING-POINT PACKAGE 
The extended floating-point package contains the following additional commands: 

0003 Floating Sine 

Take the sine of the angle in the floating accumulator (considered to be in 
radians) and leave the result in the floating accumulator. 

0004 Floating Cosine 

Take the cosine of the angle in the floating accumulator (considered to be 
in radians) and leave the result in the floating accumulator. 

0005 Floating Arc Tangent 

Take the arc tangent of the number in the floating accumulator and leave 
the result (in radians) in the floating accumulator. 

0006 Floating Exponent 

Raise e to the C(FAC) power. Leave the result in the floating accumulator. 



0007 Floating Logarithm 

Take the natural logarithm of the absolute value of the number in the floating 
accumulator. 

3-11 



Sine 

The sine routine uses the following identities: 

SIN(-X)=-SIN(X) 
SIN(X)=SIN(2ttN+F)=SIN(F) 

where 0<F<2trand N is an integer: 

SIN(tt-F)=-SIN(F) 

Using these identities, F is reduced to the range -ir/2<F<y'2. 
If Y = 2F/w so that -1<F<1, then: 
SI N(X)=(C 1 F+C 3 F 3 +C 5 F 5 +C 7 F 7 +C 9 F 9 ) 



where: 



C, = 1.57079631847 
C =-0.64596371106 
C =40.07968967928 
C 7 = -0.00467376557 
C =+0.00015148419 



The sine subroutine is valid over the range; 
-2tt.2047<X<2tt.2047 



Cosine 



The cosine routine uses the identity: 
COS(X)=SIN(tt/2-X) 

It then uses the SINE routine. 



Arc Tangent 



The arc tangent routine uses the following identities: 

ARCTAN(-X)=-ARCTAN(X) 

IF X>, 1; then ARCTAN(X)=ir/2-ARCTAN(l/X) 

FOR 0<X<1 , 
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X(A +A X 2 +A^X 4 ) 
ARCTAN(X) = - — ! — 

B +B r X 2 +B 2 X 4 

where: 

^ = +0.6402481953 
A ] =+0.4229908144 



\ 



+0.0264694361 



B =+0.6402487022 
B =+0.6363779373 
B =+0.1108328778 



The result is in the range; 

-n/2<ARCTAN(X)<ir/2 
-»<X<~ 



The logarithm routine uses the identities: 
For X<1; 

LOG(X)=LOG(2 N -F) (1<F<2) 

=NLOG(2)+LOG(F) 

8 
LOG(l+Y)= [ AY* 

C=l ' 



wh 



ere: 



A ] =+0.9999964239 

A =-0.4998741238 

A 3 = +0.3317990258 

A, =+0.2407338084 

A 5 = t0. 1676540711 

A, =-0.0953293897 
6 



Logarithm 
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A_ = +0.0360884937 
Ag =-0.0064553442 
for: 

0<X<1, the identity: 
LOG(X)=-LOG(l/X) is used. 



Floating Exponent 



The exponent routine uses the identity for XX), 
e X =2 X.|og 2 e =2 N+F =2 N 2 F 

where N is an integer and <F<1 



2 F 


=1 + 2F 




2 C 

A-F+BF - 

2 

D+F 


where: 




A = 


+9.95459578 


B = 


+0.03465735903 


C = 


+617.97226053 


D = 


+87.417497202 


Log 


2 e=l. 4426950409 


If X<0 


X 

e = 


1 

-x 
e 



Example, input A, B, and output: 

Y = LOG(COS(A/B)+/A 7 B) 

SQROOT=0002 /DEFINITIONS TO ASSEMBLER 

COS=0004 

LOG=0007 

*5 

7400 

7200 

5600 
/USES EXTENDED INTERPRETER 
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*200 






KCC 






TLS 






BEGIN, 


JMS 1 5 


/I NPUT A 




JMS 1 7 


/ENTER INTERPRETER 




FPUT A 


/STORE A 




FEXT 


/EXIT 




JMS 1 5 


/INPUT B 




JMS 1 7 


/ENTER INTERPRETER 




FPUT B 


/STORE B 




FMPY A 


/A-B 




SQROOT 


/EXTRACT ROOT 




FPUT TEMP 


/STORE IT 




FGET A 


/LOAD FAC WITH A 




FDIV B 


/DIVIDE BYB 




COS 


/TAKE COSINE 




FADD TEMP 


/ADD 




LOG 


/TAKE LOG. 




FEXT 






JMS 1 6 


/OUTPUT ANSWER 




JMP BEGIN 




A, 









B, 









TEMP, 










Output Controller 
There Is an additional routine available that formats floating-point output. It is available with both the 
basic floating-point package and the extended package. 

The controller routine requires two parameters for output formatting. It is called by an effective JMS 7200. 
At this point: 

C(62) = total number of digits to be outputted. If C(62)=0, output in E format. 

C(AC) = number of digits to the right of the decimal point. If it equals 0, do 
not type a " . " . 

If the number in the floating accumulator is larger than the field width allows, "X's" will be typed. The 
sign is typed and leading 0's are suppressed. 

Example: If the contents of the floating accumulator are 678.234 (decimal). 
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C(62) C(AC) Output 



3 





+ 678 


4 





+ 678 


4 


1 


+ 678.2 


5 


2 


+ 678.23 


6 


2 


+ 678.23 


3 


2 


+ 678. 


2 





+ XX 





1 


+ 0.678234E+03 



FLOATING POINT PACKAGE VERSIONS 

Four versions of the floating-point package are available: 

DEC-08-YQ1A-PB 

This is the basic floating-point package. It consists of the input/output package and the basic arithmetic 
instructions. Its core limits are: 

5-7; 40-61; 5600-7577 

DEC-08-YQ2A-PB 

This is the basic floating-point package with the output controller, its limits are: 
7; 40-62; 5400-7577 

The output controller does not type a carriage-return/line-feed after output and, hence, the programmer 
must provide his own routine to do so. Because of the way in which the typeout routine in the floating- 
point packages are constructed, the user must construct similar typeout routines in order to avoid timing 
errors in the teleprinter. In other words, the floating-point package uses a routine similar to: 

TSF 

JMP .-1 
TLS 

to do its typing. The user must use a similar routine, i.e., one that waits for the teleprinter flag to be 
set before executing the TLS instruction. 

DEC-08-YQ3A-PB 

This is the basic floating-point package with the extended functions. In addition, two interpretive com- 
mands are provided for input and output. The additional interpretive commands are: 



3-16 



SINE 


3 


COSINE 


4 


ARCTANGENT 


5 


EXPONENTIAL 
LOG (LN) 
INPUT 1 


6 

7 

13 


OUTPUT 


14 



Its core limits are: 

7; 40-61; 4757-7577 

Using this package, the routine on page 3-14 and 3-15 for calculating: 

Y=LOG(COS(A/B)+ Za-B) may be rewritten as: 

SQROOT=2 

COS^ 

LOG=7 

INPUT=13 

OUTPUT=14 

*7 
5600 

*200 
KCC 
TLS 

BEGIN, JMS I 7 
INPUT 



/CALL INTERPRETER 
/READ A 
/STORE IT 
/READ B 



A, 
B, 
TEMP, 

$ 



FPUT A 

INPUT 

FPUT B 

FMPY A /A- B 

SQROOT 

FPUT TEMP /(A-B)**.5 

FGET A 

FDIV B 

FCOS 

FADD TEMP 

FLOG 

OUTPUT 

FEXT 

JMP BEGIN 
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DEC-08-YQ4A-PB 

This is the basic package, the output controller, and the extended functions. This package does not 
have the interpretive I/O commands of input and output. Its core limits are: 

7; 40-62; 4557-7577 
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CHAPTER 4 
PROGRAM LISTINGS 

/FLOATING POINT ARITHMETIC INTERPRETER 







*40 






0040 


0000 


EXl, 





/OPERAND STORAGE 


0041 


0000 


AC1H, 







0042 


0000 


AC1L, 







0043 


0000 


OVER1, 







0044 


0000 


EXP, 





/F.A. 


0045 


0000 


HORD, 







0046 


0000 


LORD, 







0047 


0000 


OVER2, 







0050 


0000 


EXP1, 







0051 


0000 


QUOL, 







0052 


0000 


FPAC1, 







0053 


0000 









0054 


0000 


*61 







0061 


0000 


FLAG, 
*5600 





/ARITHMETIC ERROR FLAG 


5600 


0000 


FPNT, 







5601 


7300 




CLA CLL 




5602 


3043 




DCA OVER1 




5603 


3047 




DCA OVER2 




5604 


1600 




TAD 1 FPNT 


/GET NEXT INSTRUCTION 


5605 


3253 




DCA JUMP 




5606 


1253 




TAD JUMP 




5607 


0263 




AND PAGE NO 


/GET PAGE BIT 


5610 


7650 




SNA CLA 


/PAGE ZERO? 


5611 


5214 




JMP .+3 


Aes 


5612 


1261 




TAD MASK5 


/NO 


5613 


0200 




AND FPNT 


/C(FPNT)0-4 CONTAINS PAGE BITS 


5614 


3256 




DCA ADDR 




5615 


1262 




TAD MASK7 


/GET 7 BIT ADDRESS 


5616 


0253 




AND JUMP 




5617 


1256 




TAD ADDR 




5620 


3256 




DCA ADDR 




5621 


1264 




TAD INDRCT 


/INDIRECT BIT=1? 


5622 


0253 




AND JUMP 




5623 


7650 




SNA CLA 




5624 


5227 




JMP LOOP01 


/NO-GO ON 


5625 


1656 




TAD 1 ADDR 


/YES DEFER 


5626 


3256 




DCA ADDR 




5627 


2200 


LOOP01 , 


ISZ FPNT 




5630 


1656 




TAD 1 ADDR 




5631 


3040 




DAC EXl 


/EXPONENT 


5632 


1256 




TAD ADDR 
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5633 

5634 

5635 

5636 

5637 

5640 

5641 

5642 

5643 

5644 

5645 

5646 

5647 

5650 

5651 

5652 

5653 

5654 

5655 

5656 

5657 

5660 

5661 

5662 

5663 

5664 

5665 

5666 

5667 

5670 

5671 

5672 

5673 

5674 

5675 

5676 

5677 

5700 

5701 

5702 

5703 

5704 

5705 

5706 

5707 

5710 

5711 

5712 

5713 

5714 

5715 

5716 

5717 

5720 



3257 

2257 

1657 

3041 

2257 

1657 

3042 

1253 

7106 

7006 

0260 

1265 

3254 

1654 

3254 

5654 

0000 

0000 

0000 

0000 

0000 

0017 

7600 

0177 

0200 

0400 

5666 

5742 

5716 

5737 

5761 

6305 

5676 

5705 

5773 

1040 

3044 

1041 

3045 

1042 

3046 

5201 

1044 

3656 

2256 

1045 

3656 

2256 

1046 

3656 

5201 

4771 

5201 

4772 



JUMP, 

JUMP2, 

G02, 

ADDR, 

SAVE, 

MASK3, 

MASK5, 

MASK7, 

PAGENO, 

INDRCT, 

TABLE, 



FLGT, 



FLPT, 



FLAD, 



DCA SAVE 

ISZ SAVE 

TAD I SAVE 

DCA AC1H 

ISZ SAVE 

TAD I SAVE 

DAC AC1L 

TAD JUMP 

CLL RTL 

RTL 

AND MASK3 

TAD TABLE 

DCA JUMP2 

TAD I JUMP2 

DCA JUMP2 

JMP I JUMP2 











0017 

7600 

0177 

0200 

0400 

.+1 

EXIT 

FLAD 

FLSU 

FLMY 

FLDV 

FLGT 

FLPT 

NORF 

TAD EX1 

DCA EXP 

TAD AC1H 

DCA HORD 

TAD AC1L 

DCA LORD 

JMP FPNT+1 

TAD EXP 

DCA I ADDR 

ISZ ADDR 

TAD HORD 

DCA I ADDR 

ISZ ADDR 

TAD LORD 

DCA I ADDR 

JMP FPNT+1 

JMS I ALGN 

JMP FPNT+1 

JMS I UNORM 



/HIGH ORDER MANTISSA 



/LOW ORDER MANTISSA 



/GET BITS 0-2, IE OPCODE 
/LOOKUP IN TABLE 



/GO THERE 



/TABLE USED IN INTERPRETING 
/BITS 0-2 OF PSEUDO 
/INSTRUCTION 

/IF OPCODE=0, GO TO EXIT 
/AND INTERPRET BITS 8-11 



/FGET=5 



/FPUT=6 



/FLAD=1 - FIRST ALIGN EXPONENTS 
/RETURN IF NO ALIGNMENT IS POSSIBLE 
/LARGER OF THE TWO IS IN F.A. 
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5721 


7100 




CLL 




5722 


1043 




TAD OVER1 


/TRIPLE PRECISION ADDITION 


5723 


1047 




TAD OVER2 


/SINCE BITS ARE SHIFTED 


5724 


3047 




DCA OVER2 


/RIGHT 


5725 


7004 




RAL 




5726 


1042 




TADAC1L 




5727 


1046 




TAD LORD 




5730 


3046 




DCA LORD 




5731 


7004 




RAL 




5732 


1041 




TADAC1H 




5733 


1045 




TADHORD 




5734 


3045 




DCA HORD 




5735 


4770 




JMS 1 NORM 


/NORMALIZE THE RESULT 


5736 


5201 




JMP FPNT+1 




5737 


4741 


FLSU, 


JMS 1 OPMINS 


/FSUB=2 - NEGATE THE OPERAND 


5740 


5316 




JMP FLAD 


/ADD 


5741 


6400 


OPMINS, 


MINUS2 




5742 


1253 


EXIT, 


TAD JUMP 


/OPCODE=0 


5743 


0260 




AND MASK3 


/ARE BITS 8-11=0? 


5744 


7450 




SNA 




5745 


5600 




JMP 1 FPNT 


/YES=FEXT 


5746 


1360 




TAD ACON6 


/LOOKUP ON TABLE 


5747 


3254 




DCA JUMP2 




5750 


1654 


/ 


TAD 1 JUMP2 




5751 


3254 




DCA JUMP2 




5752 


1200 




TAD FPNT 




5753 


3255 




DCA G02 




5754 


4654 




JMS 1 JUMP2 


/CALL AS A SUBROUTINE 


5755 


1255 




TAD G02 


/RESTORE F.P. POINTER 


5756 


3200 




DCA FPNT 




5757 


5201 




JMP FPNT+1 


/GET NEXT PSEUDO INSTRUCTS 


5760 


6544 


ACON6, 


TABLE6-1 


/CALLING CAN BE TO A DEPTH O 


5761 


7201 


FLMY, 


CLA IAC 


/FMPY=3 


5762 


1040 




TAD EX1 




5763 


1044 




TAD EXP 


/ADD EXPONENTS TOGETHER 


5764 


3044 




DCA EXP 




5765 


4767 




JMS 1 MULT 


/MULTIPLY 


5766 


5201 




JMP FPNT+1 




5767 


6221 


MULT, 


DMULT 




5770 


6600 


NORM, 


DNORM 




5771 


6020 


ALGN, 


ALIGN 




5772 


6564 


UNORM, 


DUNORM 




5773 


4770 


NORF, 


JMS 1 NORM 




5774 


5201 


*6000 


JMP FPNT+1 




6000 


0000 


ACMINS, 





/ROUTINE TO PERFORM 


6001 


7300 




CLL CLA 




6002 


1047 




TAD OVER2 


/TRIPLE PRECISION NEGATION 


6003 


7041 




CMA IAC 


/OF FLOATING AC 


6004 


3047 




DCA OVER2 




6005 


1046 




TAD LORD 




6006 


7040 




CMA 




6007 


7430 




SZL 
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6010 


7101 


6011 


3046 


6012 


1045 


6013 


7040 


6014 


7430 


6015 


7101 


6016 


3045 


6017 


5600 


6020 


0000 


6021 


1045 


6022 


7640 


6023 


5227 


6024 


1046 


6025 


7650 


6026 


5337 


6027 


1041 


6030 


7640 


6031 


5235 


6032 


1042 


6033 


7650 


6034 


5620 


6035 


1040 


6036 


7041 


6037 


1044 


6040 


7450 


6041 


5314 


6042 


7500 


6043 


7041 


6044 


3342 


6045 


1342 


6046 


1343 


6047 


7700 


6050 


5256 


6051 


4316 


6052 


7430 


6053 


1350 


6054 


1347 


6055 


5325 


6056 


4316 


6057 


7420 


6060 


1350 


6061 


1347 


6062 


3344 


6063 


1342 


6064 


7041 


6065 


1744 


6066 


3744 


6067 


2344 


6070 


1344 


6071 


3345 


6072 


2345 


6073 


1345 


6074 


3346 


6075 


2346 



ALIGN, 



CLL I AC 
DCA LORD 
TAD HORD 
CMA 
SZL 

CLL I AC 
DCA HORD 

JMP I ACMINS 


TAD HORD 
SZA CLA 
JMP .+4 
TAD LORD 
SNA CLA 
JMP NOHERE 
TAD AC1H 
SZA CLA 
JMP .+4 
TADAC1L 
SNA CLA 
JMP I ALIGN 
TAD EX1 
CMA IAC 
TAD EXP 
SNA 

JMP DONE 
SMA 

CMA IAC 
DCA AMOUNT 
TAD AMOUNT 
TAD TEST2 
SMA CLA 
JMP .+6 
JMS OUTGO 
SZL 

TAD TAG2 
TAD TAG1 
JMP NOGO 
JMS OUTGO 
SNL 

TAD TAG2 
TAD TAG1 
DCA TEST3 
TAD AMOUNT 
CMA IAC 
TAD I TEST3 
DCA I TEST3 
ISZ TEST3 
TAD TEST3 
DCA TEST4 
ISZ TEST4 
TAD TEST4 
DCA TEST5 
ISZ TEST5 



/subroutine to align 
/binary points for 
/add-subtract 

/is mantissa zero? 

Aes 

/is operand zero? 



/BOTH ARE ZERO-EXIT 



/are exponents equal? 
Aes 



/can exponents be aligned? 
Aes 

/NO 



/SET UP ADDRESSES 
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6076 

6077 

6100 

6101 

6102 

6103 

6104 

6105 

6106 

6107 

6110 

6111 

6112 

6113 

6114 

6115 

6116 

6117 

6120 

6121 

6122 

6123 

6124 

6125 

6126 

6127 

6130 

6131 

6132 

6133 

6134 

6135 

6136 

6137 

6140 

6141 

6142 

6143 

6144 

6145 

6146 

6147 

6150 

6151 

6152 

6153 

6154 

6155 

6156 

6157 

6160 

6161 

6162 



7100 

1744 

7510 

7020 

7010 

3744 

1745 

7010 

3745 

1746 

7010 

3746 

2342 

5276 

2220 

5620 

0000 

1040 

7041 

1044 

7004 

7200 

5716 

3344 

1744 

3044 

2344 

1744 

3045 

2344 

1744 

3046 

5620 

1040 

3044 

5314 

0000 

0030 

0000 

0000 

0000 

0044 

7774 

5601 

1362 

3044 

1362 

3045 

7040 

3046 

2061 

5751 

3777 



SHIFT, 



DONE, 
OUTGO, 



NOGO, 



NOHERE, 



AMOUNT, 

TEST2, 

TEST3, 

TEST4, 

TEST5, 

TAG1, 

TAG2, 

RETN2, 

ERROR1, 



GOOF, 



CLL 

TAD I TEST3 

SPA 

CML 

RAR 

DCA I TEST3 

TAD I TEST4 

RAR 

DCA I TEST4 

TAD I TEST5 

RAR 

DCA I TEST5 

ISZ AMOUNT 

JMP SHIFT 

ISZ ALIGN 

JMP I ALIGN 



TAD EX1 

CMA IAC 

TAD EXP 

RAL 

CLA 

JMP I OUTGO 

DCA TEST3 

TAD I TEST3 

DCA EXP 

ISZ TEST3 

TAD I TEST3 

DCAHORD 

ISZ TEST3 

TAD I TEST3 

DCA LORD 

JMP I ALIGN 

TAD EX1 

DCA EXP 

JMP DONE 



0030 







EXP 

EX1-EXP 

FPNT+1 

TAD GOOF 

DCA EXP 

TAD GOOF 

DCA HORD 

CMA 

DCA LORD 

ISZ FLAG 

JMP I RETN2 

3777 



/THIS ROUTINE DOES 
/THE ACUTAL SHIFTING 



/DETERMINE WHICH TO SHIFT 



/CAN'T BE ALIGNED 
/LARGEST GOES INTO FAC 



/MANTISSA=0 



/DIVISION BY ZERO 
/SET TO LARGEST + VALUE 



/SET FLAG 
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6163 


0000 SQUARE, 


6164 


4407 


JMS 1 0007 


6165 


6052 


FPUT FPAC1 


6166 


3052 


FMPY FPAC1 


6167 


0000 


FEXT 


6170 


5763 


JMP 1 SQUARE 


6171 


0000 EXIT6, 


6172 


5771 

*620C 


JMP 1 EXIT6 
) 


6200 


0000 DIV1, 


o 


6201 


7300 


CLA CLL 


6202 


1045 


TAD HORD 


6203 


7510 


SPA 


6204 


7120 


CLL CML 


6205 


7010 


RAR 


6206 


3045 


DCA HORD 


6207 


1046 


TAD LORD 


6210 


7010 


RAR 


6211 


3046 


DCA LORD 


6212 


1047 


TAD OVER2 


6213 


7010 


RAR 


6214 


3047 


DCA OVER2 


6215 


7100 


CLL 


6216 


2044 


ISZ EXP 


6217 


7000 


NOP 


6220 


5600 


JMP 1 DIV1 


6221 


0000 DMULT, 


6222 


7300 


CLA CLL 


6223 


1363 


TAD SMACLA 


6224 


3347 


DCA SNSWIT 


6225 


4336 


JMS SIGN 


6226 


1042 


TAD AC1L 


6227 


3756 


DCA 1 MP2PT 


6230 


1046 


TAD LORD 


6231 


4755 


JMS 1 MP4PT 


6232 


7200 


CLA 


6233 


1757 


TAD 1 MP5PT 


6234 


3047 


DCA OVER2 


6235 


1045 


TAD HORD 


6236 


3756 


DCA 1 MP2PT 


6237 


1042 


TAD AC1L 


6240 


4755 


JMS 1 MP4PT 


6241 


1047 


TAD OVER2 


6242 


3047 


DCA OVER2 


6243 


7004 


RAL 


6244 


1757 


TAD 1 MP5PT 


6245 


3367 


DCA D 


6246 


7004 


RAL 


6247 


3370 


DCA KEEP 


6250 


1041 


TAD AC1H 


6251 


3756 


DCA 1 MP2PT 


6252 


1046 


TAD LORD 


6253 


4755 


JMS 1 MP4PT 



/DUMMY SUBROUTINE 



/SHIFT FAC RIGHT 



/DOUBLE PRECISION MULTIPLY 
/SAVE PRODUCT TRIPLE PRECISION 

/CALLS A SINGLE PRECISION 
/MULTIPLY 3 TIMES 
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6254 


1047 


6255 


3047 


6256 


7004 


6257 


1757 


6260 


1367 


6261 


3367 


6262 


7004 


6263 


1370 


6264 


3370 


6265 


1045 


6266 


3756 


6267 


1041 


6270 


4755 


6271 


1367 


6272 


3046 


6273 


7004 


6274 


1757 


6275 


1370 


6276 


3045 


6277 


4760 


6300 


3047 


6301 


2365 


6302 


5621 


6303 


4773 


6304 


5621 


6305 


1041 


6306 


7640 


6307 


5313 


6310 


1042 


6311 


7650 


6312 


5774 


6313 


1040 


6314 


7041 


6315 


1044 


6316 


7001 


6317 


3044 


6320 


1362 


6321 


3347 


6322 


4336 


6323 


4761 


6324 


1757 


6325 


1041 


6326 


7630 


6327 


7001 


6330 


3042 


6331 


3041 


6332 


2365 


6333 


4773 


6334 


5735 


6335 


5720 


6336 


0000 


6337 


1366 


6340 


3365 



TAD OVER2 
DCA OVER2 
RAL 

TAD I MP5PT 
TAD D 
DCA D 
RAL 

TAD KEEP 
DCA KEEP 
TAD HORD 
DCA I MP2PT 
TAD AC1H 
JMS I MP4PT 
TAD D 
DCA LORD 
RAL 

TAD I MP5PT 
TAD KEEP 
DCA HORD 
JMS I NORMF 
DCA OVER2 
ISZ SGN 
JMP I DMULT 
JMS I MINS 
JMP I DMULT 

FLDV, TADAC1H 

SZA CLA 
JMP .+4 
TAD AC1L 
SNA CLA 
JMP I ERROR 
TAD EX1 
CMA IAC 
TAD EXP 
IAC 

DCA EXP 
TAD SPACLA 
DCA SNSW1T 
JMS SIGN 
JMS I DIVIDE 
TAD I MP5PT 
TAD AC1H 
SZL CLA 
IAC 

DCA AC1L 
DCA AC1H 
ISZ SGN 
JMS I MINS 
JMP I .+1 
FLAD+2 

SIGN, 

TAD REST 
DCA SGN 



/DIVISION BY ZERO 



/SUBTRACT EXPONENTS 



/SET UP SIGNS 
/DIVIDE 

/ROUND OFF IN 23RD BIT 



/TEST SIGN 

/NEGATE 

/ADD ROUNDING 

/TEST SIGN OF RESULT 
/SET UP BY MULTIPLY AND 
/DIVIDE 
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6341 


1045 




TAD HORD 




6342 


7700 




SMA CLA 




6343 


5346 




JMP .+3 




6344 


4773 




JMS 1 MINS 




6345 


2365 




ISZ SGN 




6346 


1041 




TADAC1H 




6347 


7700 


SNSWIT, 


SMA CLA 


/OR SPA CLA 


6350 


5736 




JMP 1 SIGN 




6351 


4771 




JMS 1 MINS2 




6352 


2365 




ISZ SGN 




6353 


7000 




NOP 




6354 


5736 




JMP 1 SIGN 




6355 


6437 


MP4PT, 


MP4 




6356 


6471 


MP2PT, 


MP2 




6357 


6465 


MP5PT, 


MP5 




6360 


6600 


NORMF, 


DNORM 




6361 


6472 


DIVIDE, 


DUBDIV 




6362 


7710 


SPACLA, 


SPA CLA 




6363 


7700 


SMACLA, 


SMA CLA 




6364 


5601 


RETURN, 


FPNT+1 




6365 


0000 


SGN, 







6366 


7776 


REST, 


7776 




6367 


0000 


D, 







6370 


0000 


KEEP, 







6371 


6400 


MINS2, 


MINUS2 


/-AC1H,AC1L 


6372 


6420 


RAR2, 


DIV2 


/AClH,AClL/2 


6373 


6000 


MINS, 


ACM INS 




6374 


6152 


ERROR, 
*6400 


ERROR1 




6400 


0000 


MINUS2, 





/NEGATE OPERAND 


6401 


7300 




CLA CLL 


/TRIPLE PRECISION 


6402 


1043 




TAD OVER1 




6403 


7041 




CMA IAC 




6404 


3043 




DC A OVER1 




6405 


1042 




TAD AC1 L 




6406 


7040 




CMA 




6407 


7430 




SZL 




6410 


7101 




CLL IAC 




6411 


3042 




DC A AC1L 




6412 


1041 




TAD AC1H 




6413 


7040 




CMA 




6414 


7430 




SZL 




6415 


7101 




CLL IAC 




6416 


3041 




DCAAC1H 




6417 


5600 




JMP 1 MINUS2 




6420 


0000 


DIV2, 





/SHIFT OPERAND RIGHT 


6421 


7300 




CLA CLL 


ARIPLE PRECISION 


6422 


1041 




TAD AC1H 




6423 


7510 




SPA 




6424 


7120 




CLL CML 




6425 


7010 




RAR 




6426 


3041 




DCAAC1H 




6427 


1042 




TAD AC1L 





4-8 



6430 


7010 




RAR 


6431 


3042 




DC A AC1 L 


6432 


1043 




TAD OVER1 


6433 


7010 




RAR 


6434 


3043 




DCA OVER1 


6435 


7100 




CLL 


6436 


5620 




JMP 1 DIV2 


6437 


0000 


MP4, 





6440 


3266 




DCA MP1 


6441 


3265 




DCA MP5 


6442 


1270 




TAD THIR 


6443 


3267 




DCA MP3 


6444 


7100 




CLL 


6445 


1266 




TAD MP1 


6446 


7010 




RAR 


6447 


3266 




DCA MP1 


6450 


1265 




TAD MP5 


6451 


7420 




SNL 


6452 


5255 




JMP .+3 


6453 


7100 




CLL 


6454 


1271 




TAD MP2 


6455 


7010 




RAR 


6456 


3265 




DCA MP5 


6457 


2267 




ISZ MP3 


6460 


5245 




JMP MP4+6 


6461 


1266 




TAD MP1 


6462 


7010 




RAR 


6463 


7100 




CLL 


6464 


5637 




JMP 1 MP4 


6465 


0000 


MP5, 





6466 


0000 


MP1, 





6467 


0000 


MP3, 





6470 


7764 


THIR, 


7764 


6471 


0000 


MP2, 





6472 


0000 


DUBDIV, 





6473 


7300 




CLA CLL 


6474 


3051 




DCA QUOL 


6475 


1344 




TAD MIF 


6476 


3267 




DCA MP3 


6477 


5306 




JMP DVX 


6500 


1046 


DV3, 


TAD LORD 


6501 


7004 




RAL 


6502 


3046 




DCA LORD 


6503 


1045 




TAD HORD 


6504 


7004 




RAL 


6505 


3045 




DCA HORD 


6506 


1042 


DVX, 


TAD AC1 L 


6507 


1046 




TAD LORD 


6510 


3271 




DCA MP2 


6511 


7004 




RAL 


6512 


1045 




TAD HORD 


6513 


1041 




TAD AC1H 


6514 


7420 




SNL 



/SINGLE PRECISION MULTIPLY 
/12 BITS BY 12 BITS 



/DOUBLE PRECISION DIVIDE 
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6515 


5321 


JMP DV2-1 


6516 


3045 


DCA HORD 


6517 


1271 


TAD MP2 


6520 


3046 


DCA LORD 


6521 


7200 


CLA 


6522 


1 051 DV2, 


TAD QUOL 


6523 


7004 


RAL 


6524 


3051 


DCA QUOL 


6525 


1047 


TAD OVER2 


6526 


7004 


RAL 


6527 


3047 


DCA OVER2 


6530 


2267 


1SZ MP3 


6531 


5300 


JMP DV3 


6532 


1051 


TAD QUOL 


6533 


3046 


DCA LORD 


6534 


1045 


TAD HORD 


6535 


7106 


CLL RTL 


6536 


3265 


DCA MP5 


6537 


1047 


TAD OVER2 


6540 


3045 


DCA HORD 


6541 


3047 


DCA OVER2 


6542 


1265 


TAD MP5 


6543 


5672 


JMP 1 DUBDIV 


6544 


7751 MIF, 


7751 


6545 


6163 TABLE6, SQUARE 


6546 


6656 


SQROOT 


6547 


6171 


EXIT6 


6550 


6171 


EXIT6 


6551 


6171 


EXIT6 


6552 


6171 


EXIT6 


6553 


6171 


EXIT6 


6554 


6171 


EXIT6 


6555 


6171 


EXIT6 


6556 


6171 


EXIT6 


6557 


6171 


EXIT6 


6560 


6171 


EXIT6 


6561 


6171 


EXIT6 


6562 


6171 


EXIT6 


6563 


6171 


EXIT6 


6564 


0000 DUNORM, 


6565 


4220 


JMS DIV2 


6566 


4772 


JMS 1 RARl 


6567 


2040 


ISZ EX1 


6570 


7000 


NOP 


6571 


5764 


JMP 1 DUNORM 


6572 


6200 RARl , 
*6600 


DIV1 

i 


6600 


0000 DNORM, 


6601 


7300 


CLA CLL 


6602 


3255 


DCA AMT1 


6603 


3254 


DCA SIGN1 


6604 


1045 


TAD HORD 


6605 


7510 


SPA 



/TABLE FOR INTERPRETATION 

/OF BITS 8-11 

/CONTAINS ABSOLUTE ADDRESSES 

/OF PROGRAMS CALLED AS 

/SUBROUTINES 

/EXIT6=A DUMMY OR NOP 



/SHIFT OPERAND RIGHT 



/SUBROUTINE TO NORMALIZE 
/FLOATING ACCUMULATOR 



/IS MANTISSA NEGATIVE? 
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6606 


2254 




ISZ SIGN1 


/YES 


6607 


7640 




SZA CLA 


/IS MANTISSA=0? 


6610 


5217 




JMP G06 


/NO 


6611 


1046 




TAD LORD 




6612 


7640 




SZA CLA 




6613 


5217 




JMP G06 




6614 


1047 




TAD OVER2 




6615 


7650 




SNA CLA 




6616 


5251 




JMP EXIT2 


/YES 


6617 


1254 


G06, 


TAD SIGN1 


/NO 


6620 


7640 




SZA CLA 


/NEGATIVE? 


6621 


4653 




JMS 1 NEG 


/YES 


6622 


1045 


LOP, 


TAD HORD 


/WILL SHIFT BE TOO FAR? 


6623 


7104 




RAL CLL 




6624 


7710 




SPA CLA 




6625 


5241 




JMP EXIT1 


/YES 


6626 


1047 




TAD OVER2 




6627 


7104 




CLL RAL 




6630 


3047 




DCA OVER2 




6631 


1046 




TAD LORD 


/NO SHIFT MANTISSA LEFT 


6632 


7004 




RAL 




6633 


3046 




DCA LORD 




6634 


1045 




TAD HORD 




6635 


7004 




RAL 




6636 


3045 




DCA HORD 




6637 


2255 




ISZ AMT1 


/COUNT NO. OF TIMES SHIFT 


6640 


5222 




JMP LOP 




6641 


1255 


EXIT1, 


TAD AMT1 


/CORRECT EXPONENT 


6642 


7041 




CMA IAC 




6643 


1044 




TAD EXP 




6644 


3044 




DCA EXP 




6645 


1254 




TAD SIGN1 


/NEGATIVE? 


6646 


7640 




SZA CLA 




6647 


4653 




JMS 1 NEG 


/YES 


6650 


5600 




JMP 1 DNORM 




6651 


3044 


EXIT2, 


DCA EXP 


/SET TO ZERO 


6652 


5600 




JMP 1 DNORM 




6653 


6000 


NEG, 


ACMINS 




6654 


0000 


SIGN1, 







6655 


0000 


AMT1, 







6656 


0000 


SQROOT, 





/FLOATING SQUARE ROOT 


6657 


3362 




DCA FLAG1 


/TAKES ROOT OF ABSOLUTE 


6660 


4407 




JMS 15007 




6661 


6052 




FPUT FPAC1 


/VALUE 


6662 


0000 




FEXT 


/NEWTON'S METHOD IS USED 


6663 


1045 




TAD HORD 




6664 


7710 




SPA CLA 




6665 
6666 


5351 
1044 




JMP SQEND1 
TAD EXP 


/NUMBER IS NEGATIVE 


6667 
6670 


7100 
7510 




CLL 
SPA 




6671 
6672 


7020 
7010 




CML 
RAR 




6673 


3356 




DCA ITER1 


/MAKE FIRST APPROXIMATION 
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6674 


7430 




SZL 




6675 


2356 




ISZ ITER1 




6676 


7000 




NOP 




6677 


1361 




TAD SQCON1 




6700 


3357 




DCA ITER1+1 




6701 


3360 




DCA 1TER1+2 




6702 


1053 




TAD FPAC1+1 




6703 


7640 




SZA CLA 




6704 


5310 




JMP CLCU 




6705 


1054 




TAD FPAC1+2 




6706 


7650 




SNA CLA 




6707 


5354 




JMP SQEND 


/NUMBER=0 


6710 


4407 


CLCU, 


JMS 1 0007 




6711 


5052 




FGET FPAC1 




6712 


4356 




FD1V ITER1 




6713 


1356 




FADD ITER1 




67U 


0000 




FEXT 




6715 


7240 




CLA CM A 




6716 


1044 




TAD EXP 




6717 


3044 




DCA EXP 




6720 


1044 




TAD EXP 




6721 


7041 




CMA 1AC 




6722 


1356 




TAD ITER1 




6723 


7640 




SZA CLA 


/ARE EXPONENTS EQUAL? 


6724 


5345 




JMP ROOTGO 


/NO 


6725 


1045 




TAD HORD 


/ARE HIGH-ORDER MANTISSAS EQUAL? 


6726 


7041 




CMA IAC 




6727 


1357 




TADfTERl+1 




6730 


7640 




SZA CLA 




6731 


5345 




JMP ROOTGO 


/NO 


6732 


1046 




TAD LORD 




6733 


7041 




CMA IAC 




6734 


1360 




TAD ITER1+2 


/DO LOW-ORDER MANTISSAS AGREE? 


6735 


7500 




SMA 




6736 


7041 




CMA IAC 


/WITHIN ONE BIT? 


6737 


7001 




IAC 




6740 


7710 




SPA CLA 




6741 


5345 




JMP ROOTGO 


/NO 


6742 


1362 




TAD FLAG1 




6743 


3061 




DCA FLAG 




6744 


5656 




JMP 1 SQROOT 




6745 


4407 


ROOTGO, 


JMS 1 0007 




6746 


6356 




FPUT ITER1 




6747 


0000 




FEXT 




6750 


5310 




JMP CLCU 




6751 


4653 


SQEND1, 


JMS 1 NEG 


/NEGATE FAC 


6752 


2362 




ISZ FLAG1 




6753 


5260 




JMP SQROOT+2 




6754 


3044 


SQEND, 


DCA EXP 




6755 


5656 




JMP 1 SQROOT 




6756 


0000 


ITER1, 







6757 


0000 









6760 


0000 









6761 


3015 


SQCON1, 


3015 




6762 


0000 


FLAG1, 
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ACMINS 


6000 


ACON6 


5760 


AC1H 


0041 


AC1L 


0042 


ADDR 


5656 


ALGN 


5771 


ALIGN 


6020 


AMOUNT 


6142 


AMT1 


6655 


CLCU 


6710 


D 


6367 


DIVIDE 


6361 


DIV1 


6200 


DIV2 


6420 


DMULT 


6221 


DNORM 


6600 


DONE 


6114 


DUBDIV 


6472 


DUNORM 


6o64 


DVX 


6506 


DV2 


6522 


DV3 


6500 


ERROR 


6374 


ERROR 1 


6152 


EXIT 


5742 


EXIT1 


6641 


EXIT2 


6651 


EXIT6 


6171 


EXP 


0044 


EXP1 


0050 


EX1 


0040 


FLAD 


5716 


FLAG 


0061 


FLAG1 


6762 


FLDV 


6305 


FLGT 


5676 


FLMY 


5761 



FLPT 


5705 


FLSU 


5737 


FPAC1 


0052 


FPNT 


5600 


GOOF 


6162 


G02 


5655 


G06 


6617 


HORD 


0045 


INDRCT 


5664 


ITER1 


6756 


JUMP 


5653 


JUMP2 


5654 


KEEP 


6370 


LOOP01 


5627 


LOP 


6622 


LORD 


0046 


MASK3 


5660 


MASK5 


5661 


MASK7 


5662 


MIF 


6544 


MINS 


6373 


MINS2 


6371 


MINUS2 


6400 


MP1 


6466 


MP2 


6471 


MP2PT 


6356 


MP3 


6467 


MP4 


6437 


MP4PT 


6355 


MP5 


6465 


MP5PT 


6357 


MULT 


5767 


NEG 


6653 


NOGO 


6125 


N OH ERE 


6137 


NORF 


5773 


NORM 


5770 



NORMF 


6360 


OPMINS 


5741 


OUTGO 


6116 


OVER1 


0043 


OVER2 


0047 


PAGENO 


5663 


QUOL 


0051 


RAR1 


6572 


RAR2 


6372 


REST 


6366 


RETN2 


6151 


RETURN 


6364 


ROOTGO 


6745 


SAVE 


5657 


SGN 


6365 


SHIFT 


6076 


SIGN 


6336 


SIGN1 


6654 


SMACLA 


6363 


SNSWIT 


6347 


SPACLA 


6362 


SQCON1 


6761 


SQEND 


6754 


SQEND1 


6751 


SQROOT 


6656 


SQUARE 


6163 


TABLE 


5665 


TABLE6 


6545 


TAG1 


6147 


TAG2 


6150 


TEST2 


6143 


TEST3 


6144 


TEST4 


6145 


TEST5 


6146 


THIR 


6470 


UNORM 


5772 
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/FLOATING POINT I/O ROUTINES 
/REQUIRES FLOATING POINT INTERPRETER 
/ENTRY AT 0007 



0007 5600 



*7 
FPNT, 



5600 



0044 
0045 
0046 



0052 
0053 
0054 
0055 
0056 
0057 
0060 



6767 
6770 
6771 
6772 
6773 
6774 
6775 
6776 
6777 



0000 
0000 
0000 



0000 
0000 
0000 
7777 
7777 
0000 
0000 



0000 
1056 
7S50 
5767 
1377 
4776 
5767 
7344 
0212 



*44 

EXPO NT, 
HORDER, 
LORDER, 

*52 
FPAC1 , 



SWIT1, 
SWIT2, 
CHAR, 
DSWIT, 

•6767 
PRCHAR, 






I fU*&^ firwA fh^^J^^f 








7777 

7777 






/IF = 0, NO CR-LF AFTER OUTPUT 
/IF = 0, NO LF AFTER CR IN INPUT 
/CONTAINS LAST CHARACTER READ 
•= IF NO CONVERSION TOOK PLACE 



OPUT, 
LFED, 





TAD SWIT2 

SNA CLA 

JMP I PRCHAR 

TAD LFED 

JMS I OPUT 

JMP I PRCHAR 

OUT 

0212 
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/DOUBLE PRECISION DECIMAL-BINARY 
/INPUT AND CONVERSION 
♦ 7000 



IX 



ASCII -5>g 



iwr\ 



7000 


0000 


DECONV, 







7001 


7200 




CLA 




7002 


3045 




DCA 


HORDER 


7003 


3046 




DCA 


LORDER 


7004 


3265 




DCA 


SIGN 


7005 


3266 




DCA 


DNUMBR 


7006 


4342 




JMS 


INPUT 


7007 


1336 




TAD 


PLUS 


7010 


7450 




SNA 




7011 


5217 




JMP 


DECON 


7012 


1335 




TAD 


MINUS 


7013 


7440 




SZA 




7014 


5220 




JMP 


.+4 


7015 


7240 




CLA 


CMA 


7016 


3265 




DCA 


SIGN 


7017 


4342 


DECON, 


JMS 


INPUT 


7020 


7200 




CLA 




7021 


1057 




TAD 


CHAR 


7022 


1337 




TAD 


MIN9 


7023 


7500 




SMA 




7024 


5600 




JMP 


I DECONV 


7025 


1340 




TAD 


PLUS12 


7026 


7510 




SPA 




7027 


5600 




JMP 


I DECONV 


7030 


3263 




DCA 


DIGIT 


7031 


1045 




TAD 


HORDER 


7032 


0341 




AND 


MASK 


7033 


7440 




SZA 




7034 


5217 




JMP 


DECON 


7035 


2060 




ISZ 


DSWIT 


7036 


2266 




ISZ 


DNUM3R 


7037 


4241 




JMS 


MULT10 


7040 


5217 




JMP 


DECON 


7041 


0000 


MULT10, 







7042 


1046 




TAD 


LORDER 


7043 


3261 




DCA 


LTEMP 


7044 


1045 




TAD 


HORDER 


7045 


3262 




DCA 


HTEMP 


7046 


3264 




DCA 


REMAIN 


7047 


4267 




JMS 


MULT2 


7050 


4267 




JMS 


MULT2 


7051 


4303 




JMS 


DUBLAD 


7052 


4267 




JMS 


MULT2 


7053 


1263 




TAD 


DIGIT 


7054 


3261 




DCA 


LTEMP 


7055 


3262 




DCA 


HTEMP 


7056 


4303 




JMS 


DUBLAD 


7057 


1264 




TAD 


REMAIN 


7060 


5641 




JMP 


I MULT10 


7061 


0000 


LTEMP, 







7062 


0000 


HTEMP, 







7063 


0000 


DIGIT, 







7064 


0000 


REMAIN, 







7065 


0000 


SIGN, 







7066 


0000 


DNUMBR, 







7067 


0000 


MULT2, 








/INITIALIZE MANTISSA 



/TEST FOR SIGN 



/IF-, SET SWITCH 

/IS IT A DIGIT 

/NO 

/NO 
/YES 

/OVERFLOW? 

/YES- IGNORE 

/INDEX NUMBER OF DIGITS 

/CONTINUE 

/ROUTINE TO MULTIPLY 
/DOUBLE PRECISION WORD 
/BY TEN (DECIMAL) 
/REMAIN=REMAINDER 



/CALL SUBROUTINE TO 
/MULTIPLY BY TWO 
/CALL DOUBLE ADD 

/ADD LAST DIGIT RECEIVED 



/EXIT WITH REMAINDER 
/IN AC 

/DOUBLE PRECISION WORD 

/STORAGE FOR DIGIT 

/=0 IF PLUS: =7777 IF MINUS 

/=NUMBER OF DIGITS 

/MULTIPLY LORDER, HORDER BY 2 
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7070 


7300 




CLA CLL 




7071 


1046 




TAD LORDER 




7072 


7004 




RAL 




7073 


3046 




DCA LORDER 




7074 


1045 




TAD HORDER 




7075 


7004 




RAL 




7076 


3045 




DCA HORDER 




7077 


1264 




TAD REMAIN 




7100 


7004 




RAL 




7101 


3264 




DCA REMAIN 




7102 


5667 




J MP I MULT2 




7103 


0000 


DUBLAD, 





/DOUBLE PRECISION 


7104 


7300 




CLA CLL 




7105 


1046 




TAD LORDER 




7106 


1261 




TAD LTEMP 




7107 


3046 




DCA LORDER 




7110 


7004 




RAL 




7111 


1045 




TAD HORDER 




7112 


1262 




TAD HTEMP 




7113 


3045 




DCA HORDER 




7114 


7004 




RAL 




7115 


1264 




TAD REMAIN 




7116 


3264 




DCA REMAIN 




7117 


5703 




JMP I DUBLAD 




7120 


0000 


MSIGN, 





/ROUTINE TO FORM 


7121 


7300 




CLA CLL 


/2'S COMPLEMENT IF 


7122 


2265 




ISZ SIGN 


/MINUS 


7123 


5720 




JMP I MSIGN 


/SIGN=0000: EXIT 


7124 


1046 




TAD LORDER 




7125 


7041 




CMA IAC 




7126 


3046 




DCA LORDER 




7127 


1045 




TAD HORDER 




7130 


7040 




CMA 




7131 


7430 




SZL 




7132 


7001 




IAC 




7133 


3045 




DCA HORDER 




7134 


5720 




JMP I MSIGN 




7135 


7776 


MINUS, 


253-255 


/TEST FOR SIGN 


7136 


7525 


PLUS, 


-253 




7137 


7506 


MIN9, 


-272 


/TEST FOR DIGIT 


7140 


0012 


PLUS12, 


272-260 




7141 


7600 


MASK, 


7600 


/TEST FOR OVERFLOW 






/INPUT A 


, CHARACTER. IF CR, 


TEST 






/INPUT SWITCH TO SEE IF LF 


SHOULD 






/BE TYPED. IF RUBOUT, RESTART INPUT 


7142 


0000 


INPUT, 





/INPUT A CHARACTER 


7143 


7200 




CLA 




7144 


6031 




KSF 




7145 


5344 




JMP .-1 




7146 


6036 




KRB 




7147 


3057 




DCA CHAR 




7150 


1057 




TAD CHAR 




7151 


4766 




JMS I OUTPUT 




7152 


1057 




TAD CHAR 




7153 


7450 




SNA 




7154 


5343 




JMP INPUT+1 


/IGNORE BLANKS 


7155 


1370 




TAD MRBOUT 




7156 


7450 




SNA 




7157 


5767 




JMP I RESTRT 


/RUBOUT-RESTART IN] 


7160 


1371 




TAD MI NCR 
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7161 


7650 




SNA 


CLA 


7162 


4765 




J MS 


I PRINT 


7163 


1057 




TAD 


CHAR 


7164 


5742 




JMP 


I INPUT 


7165 


6767 


PRINT, 


PR CHAR 


7166 


7344 


OUTPUT, 


OUT 




7167 


7401 


RESTRT, 


FLINTP+1 


7170 


7401 


MRBOUT, 


-377 


7171 


0162 


MI NCR, 


377- 


■215 






/FLOATING 


OUTPUT "£" FORM 






/USES: 


TSF 






/ 


JMP .-1 






/ 


TLS 






♦ 7200 






7200 


0000 


FLOUTP, 







7201 


4217 




JMS 


FOUTCN 


7202 


1324 




TAD 


BEXP 


7203 


3044 




DCA 


EXPONT 


7204 


1343 




TAD 


CHE 


7205 


4344 




JMS 


OUT 


7206 


4737 




JMS 


I FEXPPT 


7207 


1055 




TAD 


SWIT1 


7210 


7650 




SNA 


CLA 


7211 


5600 




JMP 


I FLOUTP 


7212 


1341 




TAD 


CARRTN /YES 


7213 


4344 




JMS 


OUT 


7214 


1342 




TAD 


LNFEED 


7215 


4344 




JMS 


OUT 


7216 


5600 




JMP 


I FLOUTP 






/THIS WHOLE SUBROUTINE MA 






/THE OUTPUT DIGITS : CHAN 


7217 


0000 


FOUTCN, 







7220 


7300 




CLA 


CLL 


7221 


1045 




TAD 


HORDER 


7222 


7710 




SPA 


CLA 


7223 


7220 




CLA 


CML 


7224 


1327 




TAD 


SPLUS 


7225 


7430 




SZL 




7226 


1330 




TAD 


SMINUS /NO 


7227 


4344 




JMS 


OUT 


7230 


4352 




JMS 


OUTDG 


7231 


1331 




TAD 


PERIOD 


7232 


4344 




JMS 


OUT 


7233 


7300 




CLA 


CLL 


7234 


1045 




TAD 


HORDER 


7235 


7700 




SMA 


CLA 


7236 


5242 




JMP 


FGOl 


7237 


7040 




CMA 




7240 


3733 




DCA 


I SNPT 


7241 


4732 




JMS 


I MSNPT 


7242 


7240 


FGOl , 


CLA 


CMA 


7243 


1044 




TAD 


EXPONT 


7244 


3044 




DCA 


EXPONT 


7245 


3324 




DCA 


BEXP 


7246 


1044 


FG02, 


TAD 


EXPONT 


7247 


7500 




SMA 




7250 


5263 




JMP 


FG03 


7251 


1326 




TAD 


FOUR 



/CR - SEE IF TO BE FOLLOWED 

/BY LF 

/EXIT ROUTINE 



/CONVERT MANTISSA AND OUTPUT 



/CONVERT EXPONENT AND OUTPUT 
/PRINT CR-LF? 



/NO- EXIT 



/EXIT 



CHANGE JMS OUTDG TO DCA I 10, 



/NUMBER>077 



ETC. 



/NO SET 
/YES 



LINK 



/OUTPUT "0* 
/OUTPUT ".* 



/NUMBER 
/NEGATE 



IS NEGATIVE 



/SUBTRACT 1 FROM BINARY EXPONENT 
/COMPENSATE AT FG04 

/INITIALIZE DECIMAL EXPONENT 
/IS -4<EXP0NENT<-1 

/TOO LARGE: MULTIPLY BY 1/10 
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7252 


7700 




SMA 


CLA 


7253 


5270 




JMP 


FG04 


7254 


4407 




JMS 


I FPNT 


7255 


3740 




FMPY I TENPT 


7256 


0000 




FEXT 


7257 


7240 




CLA 


CMA 


7260 


1324 




TAD 


BEXP 


7261 


3324 




DCA 


BEXP 


7262 


5246 




JMP 


FG02 


7263 


4407 


FG03, 


JMS 


I FPNT 


7264 


3372 




FMPY C.10 


7265 


0000 




FEXT 


7266 


2324 




ISZ 


BEXP 


7267 


5246 




JMP 


FG02 


7270 


3734 


FG04, 


DCA 


I DPT 


7271 


4736 




JMS 


I M2PT 


7272 


4735 




JMS 


I M10PT 


7273 


7410 




SKP 




7274 


4357 


FG05A, 


JMS 


DIVTWO 


7275 


2044 




ISZ 


EXPO NT 


7276 


5274 




JMP 


FG05A 


7277 


7450 




SNA 




7300 


5311 




JMP 


FG07 


7301 


4352 


FG06, 


JMS 


OUTDG 


7302 


1325 




TAD 


MINUS7 


7303 


3044 




DCA 


EXPO NT 


7304 


4735 


FG06A, 


JMS 


I M10PT 


7305 


4352 




JMS 


OUTDG 


7306 


2044 




ISZ 


EXPONT 


7307 


5304 




JMP 


FG06A 


7310 


5617 




JMP 


I FOUTCN 


7311 


7240 


FG07, 


CLA 


CMA 


7312 


1324 




TAD 


BEXP 


7313 


3324 




DCA 


BEXP 


7314 


1045 




TAD 


HORDER 


7315 


7640 




SZA 


CLA 


7316 


5322 




JMP 


.+4 


7317 


1046 




TAD 


LORDER 


7320 


7650 




SNA 


CLA 


7321 


3324 




DCA 


BEXP 


7322 


7240 




CLA 


CMA 


7323 


5302 




JMP 


FG06+1 


7324 


0000 


BEXP, 







7325 


7772 


MINUS7, 


7772 




7326 


0004 


FOUR, 


0004 




7327 


0253 


SPLUS, 


253 




7330 


0002 


SMINUS, 


255- 


253 


7331 


0256 


PERIOD, 


256 




7332 


7120 


MSNPT, 


MSIGN 


7333 


7065 


SNPT, 


SIGN 




7334 


7063 


DPT, 


DIGIT 


7335 


7041 


M10PT, 


MULT10 


7336 


7067 


M2PT, 


MULT2 


7337 


7522 


FEXPPT, 


FEXC 




7340 


7504 


TENPT, 


TEN 




7341 


0215 


CARRTN, 


0215 




7342 


0212 


LNFEED, 


0212 




7343 


0305 


CHE, 


305 





/TOO SMALL-TIMES TEN 
/TEN 



/ONE TENTH 



/MULTIPLY BY TWO 
/IE.SHIFT LEFT 
/MULTIPLY BY TEN 

/COMPENSATE FOR 
/BINARY EXPONENT 

/IS FIRST DIGIT A ZERO 
/YES, IGNORE 
/MULTIPLICATIONS YIELD 
/DECIMAL DIGITS AS HIGH 
/ORDER REMAINDERS 
/IE. .672X10=6+. 72.. ETC 

/7 DIGITS OUTPUT?? 
/NO: CONTINUE 
/YES: EXIT 

/IGNORE FIRST DIGIT 
/SUBTRACT 1 FROM 
/DECIMAL EXPONENT 



/IS MANTISSA ZERO? 



/YES: EXP=0 



/CONTAINS DECIMAL EXPONENT 
/NUMBER OF DIGTS OUTPUT 



/POINTERS 
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7344 
7345 


0000 
6041 


OUT, 



TSF 




7346 


5345 




JMP 


.-1 


7347 


6046 




TLS 




7350 


7200 




CLA 




7351 


5744 




JMP 


I OUT 


7352 


0000 


OUTDG, 







7353 


1356 




TAD 


C260 


7354 


4344 




JMS 


OUT 


7355 


5752 




JMP 


I OUTDG 


7356 


0260 


C260, 


0262 


1 


7357 


0000 


DIVTWO, 







7360 


7110 




CLL 


RAR 


7361 


3344 




DCA 


OUT 


7362 


1045 




TAD 


HORDER 


7363 


7010 




RAR 




7364 


3045 




DCA 


HORDER 


7365 


1046 




TAD 


LORDER 


7366 


7010 




RAR 




7367 


3046 




DCA 


LORDER 


7370 


1344 




TAD 


OUT 


7371 


5757 




JMP 


I DIVTWO 


7372 


7775 


C.10, 


7775 




7373 


3146 




3146 




7374 


3147 




3147 


r 






/FLOATING 


POINT INPUT 






*7400 






7400 


0000 


FLINTP, 







7401 


7240 




CLA 


CMA 


7402 


3313 




DCA 


PRSW 


7403 


3060 




DCA 


DSWIT 


7404 


4716 




JMS 


I DPCVPT 


7405 


7200 




CLA 




7406 


1057 




TAD 


CHAR 


7407 


1312 




TAD 


PER 


7410 


7640 




SZA 


CLA 


7411 


5220 




JMP 


FIGOl 


7412 


1313 




TAD 


PRSW 


7413 


7650 




SNA 


CLA 


7414 


5222 




JMP 


FIG02 


7415 


3721 




DCA 


I DPN 


7416 


3313 




DCA 


PRSW 


7417 


5717 




JMP 


I DPCSPT 


7420 


1313 


FIGOl , 


TAD 


PRSW 


7421 


7650 




SNA 


CLA 


7422 


1721 


FIG02, 


TAD 


I DPN 


7423 


7041 




CMA 


IAC 


7424 


3314 




DCA 


SEXP 


7425 


4720 




JMS 


I MSGNPT 


7426 


1311 


FIG03, 


TAD 


C27 


7427 


3044 




DCA 


EXPONT 


7430 


4407 




JMS 


I FPNT 


7431 


7000 




FNOR 


7432 


6052 




FPUT FPAC1 


7433 


0000 




FEXT 


7434 


1057 




TAD 


CHAR 


7435 


1310 




TAD 


MINUSE 



/OUTPUT ONE ASCII CHARACTER 



/OUTPUT ONE DIGIT 



/DIVIDE BY TWO IE. 
/ROTATE RIGHT 
/TEMPORARY STORAGE 



/CONSTANT .10 USED IN 
/FLOATING OUTPUT-PROVIDES 



/INITIALIZE "PERIOD SWITCH' 



/7777 = NO PERIOD 



/PERIOD FOUND 

/SECOND PERIOD 

/YES, TERMINATE 

/NO - SET NUMBER OF DIGITS TO 

/SET PERIOD SWITCH TO 

/CONVERT REST OF STRING 

/PERIOD READ IN PREVIOUSLY? 

/YES: -NUMBER OF DIGITS IN SEXP 
/NO 

/TEST SIGN 

/NORMALIZE F.P. NUMBER 
/SAVE NUMBER 
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7436 


7640 




SZA CLA 


/"E" READ IN? 


7437 


5252 




JMP ENDFI 


/NO 


7443 


4716 




JMS I DPCVPT 


/YES - CONVERT DECIMAL EXPONENT 


7441 


4720 




JMS I MSGNPT 


/TEST SIGN 


7442 


1045 




TAD HORDER 


/EXPONENT TOO LARGE?? 


7443 


7510 




SPA 




7444 


7001 




IAC 




7445 


7640 




SZA CLA 




7446 


5277 




JMP EXCESS 


/YES 


7447 


1046 




TAD LORDER 


/NO: DECIMAL POINT 13 


7450 


1314 




TAD SEXP 


/C( SEXP) PLACES TO RIGHT 


7451 


3314 




DCA SEXP 


/OF LAST DIGIT 






/END OF 


FLOATING POINT 


INPUT 






/COMPENSATE FOR DECIMAL 


EXPONENTS 


7452 


4407 


ENDFI, 


JMS I FPNT 


/RESTORE MANTISSA 


7453 


5052 




FGET FPAC1 




7454 


0000 




FEXT 




7455 


1314 




TAD SEXP 




7456 


7450 




SNA 




7457 


5600 




JMP I FLINTP 




7460 


7700 




SMA CLA 




7461 


5270 




JMP FIG04 




7462 


4407 




JMS I FPNT 


/. IS TO THE LEFT: 


7463 


3707 




FMPY I PC. 10 


/TIMES .1000 


7464 


0000 




FEXT 




7465 


2314 




ISZ SEXP 




7466 


5255 




JMP ENDFI+3 




7467 


5600 




JMP I FLINTP 




7470 


4407 


FIG04, 


JMS I FPNT 


/. 15 TO THE RIGHT: 


7471 


3304 




FMPY TEN 


/MULTIPLY BY 10 


7472 


0000 




FEXT 




7473 


7240 




CLA CMA 




7474 


1314 




TAD SEXP 




7475 


3314 




DCA SEXP 




7476 


5255 




JMP ENDFI+3 




7477 


1315 


EXCESS, 


TAD C3 777 




7500 


3044 




DCA EXPONT 




7501 


1315 




TAD C3777 




7502 


3045 




DCA HORDER 




7503 


5600 




JMP I FLINTP 




7504 


0004 


TEN, 


0004 




7505 


2400 




2400 




7506 


0000 




0000 




7507 


7372 


PC. 10, 


C.10 


/.10 


7510 


7473 


MINUSE, 


-305 




7511 


0027 


C27, 


0027 




7512 


7522 


PER, 


-256 




7513 


0000 


PRSW, 







7514 


0000 


SEXP, 





/CONTAINS DECIMAL EXPONENT 


7515 


3777 


C3777, 


3777 




7516 


7000 


DPCVPT, 


DECONV 




7517 


7017 


DPCSPT, 


DECON 




7520 


7120 


MSGNPT, 


MSIGN 




7521 


7066 


DPN, 
/OUTPUT 


DNUMBR 
THE EXPONENT 




7522 


0000 


FEXC, 







7523 


7300 




CLA CLL 
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7524 


1044 




TAD 


EXPONT 


7525 


7510 




SPA 




7526 


7061 




CMA 


IAC CML 


7527 


3044 




DCA 


EXPONT 


7530 


1366 




TAD 


C253 


7531 


7430 




SZL 




7532 


1367 




TAD 


C255 


7533 


4774 




J MS 


I DGPT 


7534 


3045 




DCA 


HORDER 


7535 


1044 




TAD 


EXPONT 


7536 


2045 




ISZ 


HORDER 


7537 


1370 




TAD 


Ml 44 


7540 


7500 




SMA 




7541 


5336 




JMP 


.-3 


7542 


1371 




TAD 


C144 


7543 


3044 




DCA 


EXPONT 


7544 


7040 




CMA 




7545 


1045 




TAD 


HORDER 


7546 


7440 




SZA 




7547 


4774 




JMS 


I DGPT 


7550 


3045 




DCA 


HORDER 


7551 


1044 




TAD 


EXPONT 


7552 


2045 




ISZ 


HORDER 


7553 


1372 




TAD 


Ml 2 


7554 


7500 




SMA 




7555 


5352 




JMP 


.-3 


7556 


1373 




TAD 


C12 


7557 


3046 




DCA 


LORDER 


7560 


7240 




CLA 


CMA 


7561 


1045 




TAD 


HORDER 


7562 


4774 




JMS 


I DGPT 


7563 


1046 




TAD 


LORDER 


7564 


4774 




JMS 


I DGPT 


7565 


5722 




JMP 


I FEXC 


7566 


7773 


C253, 


025^ 


>-260 


7567 


0002 


C255, 


255- 


■253 


7570 


7634 


Ml 44, 


7634 




7571 


0144 


CI 44, 


0144 


I 


7572 


7766 


M12, 


7766 




7573 


0012 


C12, 


0012 


i 


7574 


7352 


DGPT, 


0UT3G 



4-21 



BEXP 

CftRRTN 

CHAR 

CHE 

C. lii 

CI 2 

CI 44 

C253 

C255 

C260 

C27 

C5777 

DECON 

DECONV 

DGPT 

DIGIT 

DIVTWO 

DNUMSR 

DPCSPT 

DPCVPT 

DPN 

DPT 

DSWIT 

DUB LAD 

ENDFI 

EXCESS 

EXPO NT 

FEXC 

FEXPPT 

FGOl 

FG02 

FG03 

FG04 

FG05A 

FG06 

FG06A 

FG07 

FIGOl 

FIG02 

FIG 03 

FIG04 

FLINTP 

FLOUTP 

FOUR 

FOUTCN 

FPAC1 

FPNT 



7324 

7341 

0057 

7343 

73 72 

7573 

7571 

7566 

7567 

7356 

7511 

7515 

7017 

7000 

7574 

7063 

7357 

7066 

7517 

7516 

7521 

7334 

0060 

7103 

7452 

7477 

0044 

7522 

7337 

7242 

7246 

7263 

7270 

7274 

7301 

7304 

7311 

7420 

7422 

7426 

7473 

7400 

7200 

7326 

7217 

0052 

0007 



HORDER 

HTEMP 

INPUT 

LFED 

LNFEED 

LORDEK 

LTEMP 

MASK 

MI NCR 

MI NUS 

MINUSE 

MINUS7 

MIN9 

MRS OUT 

MSGNPT 

MSIGN 

MSNPT 

MULT10 

MULT2 

Ml 0PT 

Ml 2 

Ml 44 

M2PT 

OPUT 

OUT 

OUTDG 

OUTPUT 

PC 10 

PER 

PERIOD 
PLUS 

PLUS 12 

PRCHAR 

PRINT 

PRSW 

REMAIN 

RESTRT 

SEXP 

SIGN 

SMI NUS 

SNPT 

SPLUS 

SWIT1 

SWIT2 

TEN 

TENPT 



0045 

7062 

7142 

6777 

7342 

0046 

7061 

7141 

7171 

7135 

7510 

7325 

7137 

7170 

7520 

7120 

7332 

7041 

7067 

7335 

7572 

7570 

7336 

6776 

7344 

7352 

7166 

7507 

7512 

7331 

7136 

7140 

6767 

7165 

7513 

7064 

7167 

7514 

7065 

7330 

7333 

7327 

0055 

0056 

7504 

7340 
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EXTENDED FUNCTIONS 
/FLOATING POINT PACKAGE 

GETSGN=TAD 45 
CLFPNTrJMS I 7 

/FLOATING POINT EXPONENTIAL 







♦5000 




5000 


0000 


FEXP, 




5001 


1045 


GETSGN 


5002 


7710 


SPA CLA 


5003 


4310 


JMS FNEG 


5004 


3304 


DCA SIGN 


5005 


4407 


CLFPNT 


5006 


3273 


FMPY 


LG2E 


5007 


6315 


FPUT 


X 


5010 


0000 


FEXT 




5011 


4706 


JMS I FIXIT 


5012 


1045 


GETSGN 


5013 


3305 


DCA FLAG2 


5014 


4707 


JMS I FLOAT 


5015 


4407 


CLFPNT 


5016 


6320 


FPUT 


XSQR 


5017 


5315 


FGET 


X 


5020 


2320 


FSUB 


XSQR 


5W21 


6315 


FPUT 


X 


5022 


3315 


FMPY 


X 


5023 


6320 


FPUT 


XSQR 


5024 


1270 


FADD 


D 


5025 


6323 


FPUT 


TEMP 


5026 


5265 


FGET 


C 


5027 


4323 


FDIV 


TEMP 


5030 


2315 


FSUB 


X 


5031 


1257 


FADD 


A 


5032 


6323 


FPUT 


TEMP 


5033 


5262 


FGET 


B 


5034 


3320 


FMPY 


XSQR 


5035 


1323 


FADD 


TEMP 


5036 


6323 


FPUT 


TEMP 


5037 


5315 


FGET 


X 


5040 


4323 


FDIV 


TEMP 


5041 


3301 


FMPY 


TWO 


5042 


1276 


FADD 


ONE 


5043 


0000 


FEXT 




5044 


1305 


TAD FLAG2 


5045 


1044 


TAD 44 


5046 


3044 


DCA 44 


5047 


2304 


ISZ SIGN 


5050 


5600 


JMP I FEXP 


5051 


4407 


CLFPNT 


5052 


6315 


FPUT 


X 


5053 


5276 


FGET 


ONE 


5054 


4315 


FDIV 


X 


5055 


0000 


FEXT 




5056 


5600 


JMP I FEXP 






/CONSTANTS FOR 


FEXP 


5057 


0004 


A, 0004 




5060 


2372 


2372 




5061 


1402 


1402 





/C(SIGN)=-1 IF X<0 
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5062 


7774 


B, 


7774 


5063 


2157 




2157 


5064 


5157 




5157 


5065 


0012 


c, 


0012 


5066 


5454 




5454 


5067 


0343 




0343 


5070 


0007 


D, 


0007 


5071 


2566 




2566 


5072 


5341 




5341 


5073 


0001 


LG2E, 


0001 


5074 


2705 




2705 


5075 


2435 




2435 


5076 


0001 


ONE, 


0001 


5077 


2000 




2000 


5100 


0000 




0000 


5101 


0002 


TWO, 


0002 


5102 


2000 




2000 


5103 


0000 




0000 


5104 


0000 


SIGN, 





5105 


0000 


FLAG2, 





5106 


4757 


FIXIT, 


FIX 


5107 


5563 


FLOAT, 


FLOA 






/NEGATION 


SUBROUTINE 


5110 


0000 


FNEG, 





5111 


4714 




JMS I ACMINS 


5112 


7240 




CLA CMA 


5113 


5710 




JMP I FNEG 


5114 


6000 


ACMINS, 


6000 






/TEMPORARY STORAGE 


5115 


0000 


X, 





5116 


0000 







5117 


0000 







5120 


0000 


XSQR, 





5121 


0000 







5122 


0000 







5123 


0000 


TEMP, 





5124 


0000 







5125 


0000 







5126 


0000 


ZER, 





5127 


0000 







5130 


0000 











/MAIN ALGORITHM FOR ARCTANGENT 


5131 


4407 


ARCALG, 


CLFPNT 


5132 


5315 




FGET X 


5133 


3315 




FMPY X 


5134 


6320 




FPUT XSQR 


5135 


3373 




FMPY BET2 


5136 


1370 




FADD BET1 


5137 


3320 




FMPY XSQR 


5140 


1365 




FADD BETZ 


5141 


6323 




FPUT TEMP 


5142 


5362 




FGET ALF2 


5143 


3320 




FMPY XSQR 


5144 


1357 




FADD ALF1 


5145 


3320 




FMPY XSQR 


5146 


1354 




FADD ALFZ 


5147 


3315 




FMPY X 


5150 


4323 




FDIV TEMP 


5151 


0000 




FEXT 


5152 


5753 




JMP I .+1 


5153 


5226 




ARCRTN 



/CALL SUBROUTINE IN 
/INTERPRETER 

/POINTER TO INTERPRETER 
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/CONSTANTS - FLOATING ARC TANGENT 



5154 


0000 


ALFZ, 


0000 




5155 


2437 




2437 




5156 


1643 




1643 




5157 


7777 


ALF1, 


7777 




5160 


3304 




3304 




5161 


4434 




4434 




5162 


7773 


ALF2, 


7773 




5163 


3306 




3306 




5164 


5454 




5454 




5165 


0000 


BETZ, 


0000 




5166 


2437 




2437 




5167 


1646 




1646 




5170 


0000 


BET1, 


0000 




5171 


2427 




2427 




5172 


2323 




2323 




5173 


7775 


BET2, 


7775 




5174 


3427 




3427 




5175 


7052 


/FLOATING 
*5200 


7052 
i POINT ARC TANGENT 




5200 


0000 


ARTN, 







5201 


1045 




GETSGN 




5202 


7710 




SPA CLA 




5203 


4642 




JMS I NEGIT 




5204 


3243 




DCA FLAG1 




5205 


4407 




CLFPNT 




5206 


6644 




FPUT I XI 




5207 


5644 




FGET I XI 




5210 


2645 




FSUB I CON1 




5211 


0000 




FEXT 




5212 


1045 




GETSGN 




5213 


7710 




SPA CLA 




5214 


5223 




JMP GO 




5215 


4407 




CLFPNT 




5216 


5645 




FGET I CON1 




5217 


4644 




FDIV I XI 




5220 


6644 




FPUT I XI 




5221 


0000 




FEXT 




5222 


7240 




CLA CMA 




5223 


3251 


GO, 


DCA FLOG 


/FLAG 


5224 


5625 




JMP I .+1 


/CALL ALGORITHM 


5225 


5131 




ARCALG 




5226 


2251 


ARCRTN, 


ISZ FLOG 


/RETURN HERE 


5227 


5235 




JMP EXIT 




5230 


4407 




CLFPNT 




5231 


6644 




FPUT I XI 




5232 


5246 




FGET PIOT 




5233 


2644 




FSUB I XI 




5234 


0000 




FEXT 




5235 


2243 


EXIT, 


ISZ FLAG1 




5236 


5600 




JMP I ARTN 




5237 


4642 




JMS I NEGIT 




5240 


7200 




CLA 




5241 


5600 




JMP I ARTN 








/CONSTANTS FOR ARCTANGENT 




5242 


5110 


NEGIT, 


FNEG 




5243 


0000 


FLAG1, 







5244 


5115 


XI, 


X 




5245 


5076 


CON1, 


ONE 
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5246 


0001 


PIOT, 0001 






5247 


3110 


3110 






5250 


3755 


3755 










/FLOATING LOGARITHM 




5251 


0000 


FLOG, 






5252 


1045 


GETSGN 




5253 


7440 


SZA 






5254 


5261 


JMP OK 




5255 


4407 


CLFPNT 


/INDICATE ERROR 


5256 


4775 


FDIV 


I ZERO 


/DIVIDE BY ZERO 


5257 


0000 


FEXT 






5260 


5651 


JMP I FLOG 




5261 


7710 


OK, SPA CLA 




5262 


4642 


JMS ] 


: NEGIT 




5263 


4407 


CLFPNT 




5264 


6774 


FPUT 


I TEM 




5265 


2645 


FSUB 


I CON1 




5266 


0000 


FEXT 






5267 


1045 


GETSGN 




5270 


7450 


SNA 






5271 


5350 


JMP ZERGO 




5272 


7710 


SPA CLA 




5273 


5341 


JMP INVERT 




5274 


3243 


START, DCA FLAG1 




5275 


7040 


CMA 






5276 


1774 


TAD ] 


[ TEM 




5277 


3045 


DCA 45 




5300 


4777 


JMS I 


: FLOATP 




5501 


4407 


CLFPNT 




5502 


3776 


FMPY 


I L0G2 




5305 


6644 


FPUT 


I XI 




55u4 


0000 


FEXT 






55id5 


7001 


IAC 






5306 


3774 


DCA I 


. TEM 




5307 


4407 


CLFPNT 




5310 


5774 


FGET 


I TEM 




5311 


2645 


FSUB 


I CON1 




5312 


6774 


FPUT 


I TEM 




5313 


3773 


FMPY 


I L8 




5314 


1772 


FADD 


I L7 




5315 


3774 


FMPY 


I TEM 




5316 


1771 


FADD 


I L6 




5317 


3774 


FMPY 


I TEM 




5320 


1770 


FADD 


I L5 




5321 


3774 


FMPY 


I TEM 




5322 


1365 


FADD 


L4 




5323 


3774 


FMPY 


I TEM 




5324 


1362 


FADD 


L3 




5325 


3774 


FMPY 


I TEM 




5326 


1357 


FADD 


L2 




5327 


3774 


FMPY 


I TEM 




5330 


1354 


FADD 


LI 




5331 


3774 


FMPY 


I TEM 




5332 


1644 


FADD 


I XI 




5333 


0000 


FEXT 






5334 


2243 


ISZ FLAG1 




5335 


5651 


JMP ] 


: FLOG 




5336 


4642 


JMS ] 


[ NEGIT 




5337 


7200 


CLA 






5340 


5651 


JMP ] 


[ FLOG 
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5341 


4407 


INVERT, 


CLFPNT 




5342 


5645 




FGET I CON1 




5343 


4774 




FDIV I TEM 




5344 


6774 




FPUT I TEM 




5345 


0000 




FEXT 




5346 


7240 




CLA CMA 




5347 


5274 




JMP START 




5350 


4407 


ZERGO, 


CLFPNT 




5351 


5775 




FGET I ZERO 




5352 


0000 




FEXT 




5353 


5651 




JMP I FLOG 




5354 


0000 


LI, 


0000 




5355 


3777 




3777 




5356 


7742 




7742 




5357 


7777 


L2, 


7777 




5360 


4000 




4000 




5361 


4100 




4100 




5362 


7777 


L3, 


7777 




5363 


2517 




2517 




5364 


0307 




0307 




5365 


7776 


L4, 


7776 




5366 


4113 




4113 




5367 


7211 




7211 




5370 


5547 


L5, 


L0G5 




5371 


5552 


L6, 


L0G6 




5372 


5555 


L7, 


LOG 7 




5373 


5560 


LS, 


LOGS 




5374 


5123 


TEM, 


TEMP 




5375 


5126 


ZERO, 


ZER 




5376 


5544 


L0G2, 


L0GE2 




5377 


5563 


FLOATP, 

/FLOATING 

♦5400 


FLOA 
POINT SINE 




5400 


0000 


FSIN, 







5401 


1045 




GETSGN 


/X>0? 


5402 


7740 




SMA SZA CLA 




5403 


5210 




JMP MOD 


/YES 


5404 


1045 




GETSGN 




5405 


7700 




SMA CLA 


/NO X=0? 


5406 


5600 




JMP I FSIN 


/YES SIN(0)=0 


5407 


4741 




JMS I NEGT 


/NO:SIN(-X)=-SIN(X) 


5410 


3343 


MOO, 


DCA PNTR 


/REDUCE X MODULO 2 PI 


5411 


4407 




CLFPNT 




5412 


4315 




FDIV TWOPI 




5413 


6724 




FPUT I XSQ2 




5414 


0000 




FEXT 




5415 


4742 




JMS I FIXR 




5416 


4363 




JMS FLOA 




5417 


4407 




CLFPNT 




5420 


6723 




FPUT I X2 




5421 


5724 




FGET I XSQ2 




5422 


2723 




FSUB I X2 




5423 


3315 




FMPY TWOPI 




5424 


6723 




FPUT I X2 




5425 


2320 




FSUB PI 


/X<PI? 


5426 


0000 




FEXT 




5427 


1045 




GETSGN 




5430 


7710 




SPA CLA 




5431 


5241 




JMP PCHECK 


/YES 


5432 


4407 




CLFPNT 


/NO SIN(X-PI)=-SIN(X> 


5433 


6723 




FPUT I X2 




5434 


0000 




FEXT 
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5435 1343 

5436 7653 

5437 7040 

5440 3343 

5441 4407 PCHECK, 

5442 5723 

5443 2714 

5444 0000 

5445 1045 

5446 7710 

5447 5255 

5450 4407 

5451 5320 

5452 2723 

5453 6723 

5454 0000 



X2 
PI2 



TAD PNTR 
SNA CLA 
CMA 

DCA PNTR 
CLFPNT 
FGET I 
FSUB I 
FEXT 
GETSGN 
SPA CLA 
JMP PALG 
CLFPNT 
FGET PI 
FSUB I X2 
FPUT I X2 
FEXT 



/X<PI/2? 



/YES 

/NO 

/SIN(X)=SIN(PI-X) 



5455 
5456 
5457 
5460 
5461 
5462 
5463 
5464 
5465 
5466 
5467 
5470 
5471 
5472 
5473 
5474 
5475 
5476 
5477 
5500 
5501 
5502 



4407 
5723 
4714 
6723 
3723 
6724 
5325 
3724 
1330 
3724 
1333 
3724 
1336 
3724 
1714 
3723 
0000 
2343 
5600 
4741 
7200 
5600 



PALG, 



CLFPNT 



FGET 
FDIV 
FPUT 
FMPY 
FPUT 



X2 

PI2 

X2 

X2 

XSQ2 



FGET C9 

FMPY I XSQ2 
FADD C7 

FMPY I XSQ2 
FADD C5 
FMPY I XSQ2 
FADD C3 



FMPY 
FADD 
FMPY 
FEXT 
ISZ PNTR 



I XSQ2 
I PI2 
I X2 



JMP 
JMS 
CLA 
JMP 



I FSIN 
I NEGT 

I FSIN 



/FLOATING POINT COSINE 



5503 
5504 
5505 
5506 
5507 
5510 
5511 
5512 
5513 



0000 
4407 
6723 
5714 
2723 
0000 
1303 
3200 
5201 



FCOS, 





CLFPNT 
FPUT I 
FGET I 
FSUB I 
FEXT 

TAD FCOS 
DCA FSIN 
JMP FSIN+1 



X2 

PI2 

X2 



/C0S(X)=SIN(PI/2-X) 



/CONSTANTS AND POINTERS 



5514 
5515 
5516 
5517 
5520 
5521 
5522 
5523 
5524 



5246 
0003 
3110 
3755 
0002 
3110 
3755 
5115 
5120 



PI2, 
TWOPI, 



PI. 



X2, 
XSQ2, 



PIOT 

0003 

3110 

3755 

0002 

3110 

3755 

X 

XSQR 



/PI/2 
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/SINE 


CONSTANTS 


5525 


7764 


C9, 


7764 


5526 


2366 




2366 


5527 


5735 




5735 


5530 


7771 


C7, 


7771 


5531 


5466 




5466 


5532 


6317 




6317 


5533 


7775 


C5, 


7775 


5534 


2431 




2431 


5535 


5053 




5053 


5536 


0000 


C3, 


0000 


5537 


5325 




5325 


5540 


0420 




0420 


5541 


5110 


NEGT, 


FNEG 


5542 


4757 


FIXR, 


FIX 


5543 


0000 


PNTR, 









/LOGARITHM CONSTANTS 


5544 


0000 


L0GE2, 


0000 


5545 


2613 




2613 


5546 


4414 




4414 


5547 


7776 


LOG5, 


7776 


5550 


2535 




2535 


5551 


3301 




3301 


5552 


7775 


L0G6, 


7775 


5553 


4746 




4746 


5554 


0771 




0771 


5555 


7774 


LOG7, 


7774 


5556 


2236 




2236 


5557 


4304 




4304 


5560 


7771 


LOG 8, 


7771 


5561 


4544 




4544 


5562 


1735 




1735 






/FLOAT 


CC45) 


5563 


0000 


FLOA, 





5564 


7300 




CLA CLL 


5565 


3046 




OCA 46 


5566 


1374 




TAD CI 3 


55S7 


3044 




DCA 44 


5570 


4407 




CLFPNT 


5571 


7000 




FNOR 


5572 


0000 




FEXT 


5573 


5763 




JMP I FLOA 


5574 


0013 


C13, 
♦4757 


0013 






/FIX C(FAC) 


4757 


0000 


FIX, 





4760 


1044 




TAD 44 


4761 


7540 




SMA SZA 


4762 


5365 




JMP .+3 


4763 


7200 




CLA 


4764 


5375 




JMP FIXEND 


4765 


1377 




TAD Ml 3 


4766 


3044 




DCA 44 


4767 


1044 


LOOP, 


TAD 44 


4770 


7700 




SMA CLA 


4771 


5757 




JMP I FIX 


4772 


4774 




JMS I .+2 


4773 


5367 




JMP LOOP 


4774 


6200 




6200 



/DIV1 IN INTERPRETER 
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4775 
4776 
4777 



3045 
5757 
7765 



FIXEND, 
M13, 



DCA 45 
JMP I FIX 
-13 



A 


5057 


ACrilNS 


5114 


ALFZ 


5154 


ALF1 


5157 


ALF2 


5162 


ARCALG 


5131 


ARCRTN 


5226 


ARTN 


5200 


B 


5062 


BETZ 


5165 


BETl 


5170 


BET2 


5173 


C 


5065 


CLFPNT 


4407 


C0N1 


5245 


CI 3 


5574 


C3 


5536 


C5 


5533 


C7 


5530 


C9 


5525 


D 


5070 


EXIT 


5235 


FCOS 


5503 


FEXP 


5000 


FIX 


4757 


FIXEND 


4775 


FIXIT 


5106 


FIXR 


5542 


FLAG1 


5243 


FLAG2 


5105 


FLOA 


5563 


FLOAT 


5107 


FLOATP 


5377 


FLOG 


5251 


FNEG 


5110 


F5IN 


54d0 


GtTSGN 


1045 


GO 


5223 


INVERT 


5341 


LG2E 


5073 


L0GE2 


5544 



L0G2 


5376 


LOG 5 


5547 


L0G6 


5552 


LOG7 


5555 


LOG8 


5560 


LOOP 


4767 


LI 


5354 


L2 


5357 


L3 


5362 


L4 


5365 


L5 


5370 


L6 


5371 


L7 


5372 


L8 


5373 


MOD 


5410 


Ml 3 


4777 


NEGIT 


5242 


NEGT 


5541 


OK 


5261 


ONE 


5076 


PALG 


5455 


PCriECK 


5441 


PI 


5520 


PIOT 


5246 


PI 2 


5514 


PNTK 


5543 


SIGN 


5104 


STAKT 


5274 


TEM 


5374 


TEMP 


5123 


TWO 


5101 


TVOPI 


5515 


X 


5115 


X5QR 


5120 


XSQ2 


5524 


XI 


5244 


X2 


5523 


ZtK 


5126 


ZEKGO 


5350 


ZERO 


5375 
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/FLOATING OUTPUT PROGRAM 
/IF(62>=0, THEN OUTPUT IN E FORMAT 
/ELSE, C(62)=NUMBER OF DIGITS 
/CCAC)=NUMBER OF PLACES TO RIGHT OF . 
/IF C(AC>=0, THEN DON'T OUTPUT POINT 
/SIGN AND . NOT COUNTED IN OUTPUT 
/CONTENTS OF 15 LOST DURING OPERATION 



7201 
7202 
7203 
7204 



4777 
4217 
1324 
4776 



*7201 



JMS I 7377 /ALTERATIONS FLOATING OUTPUT 
JMS 7217 
TAD 7324 
JMS I 7376 



720 7 7200 



*7207 



CLA 



*7227 



7227 


3415 




DCA 


I 15 


7230 


7000 




NOP 




7231 


7000 




NOP 




7232 


7000 


♦ 7301 


NOP 




7301 


3415 


♦ 7305 


DCA 


I 15 


7305 


3415 


♦ 7376 


DCA 


I 15 


7376 


5412 




TGO 




7377 


5400 




EDIT 






♦ 5400 






5400 


0000 


EDIT, 







5401 


3352 




DCA 


SAC 


5402 


1062 




TAD 


62 


5403 


7041 




CMA 


IAC 


5404 


3347 




DCA 


COUNT1 


5405 


1357 




TAD 


M8 


5406 


3350 




DCA 


C0UNT2 


5407 


1363 




TAD 


SAD1 


5410 


3015 




DCA 


15 


5411 


5600 




J MP 


I EDIT 


5412 


0000 


TGO, 







5413 


3044 




DCA 


44 


5414 


1363 




TAD 


SAD1 


5415 


3015 




DCA 


15 


5416 


1415 




TAD 


I 15 


5417 


4761 




JMS 


i oun 


5420 


1062 




TAD 


62 


5421 


7650 




SNA 


CLA 


5422 


5327 




J MP 


EFORM 


5423 


2212 




ISZ 


TGO 


5424 


2212 




ISZ 


TGO 


5425 


1044 


TRYAGN, 


TAD 


44 


5426 


7510 




SPA 




5427 


5300 




J MP 


MINS 


5430 


1352 




TAD 


SAC 


5431 


7041 




CMA 


IAC 


5432 


1062 




TAD 


62 
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5433 


7510 


SPA 




5434 


5266 


JMP 


ERR 


5435 


7450 


SNA 




5436 


5245 


JMP 


G01 


5437 


7041 


CMA 


IAC 


5440 


3351 


DCA 


CNTR 


5441 


1353 


TAO 


SPCE 


5442 


4321 


JMS 


OUT 


5443 


2351 


ISZ 


CNTR 


5444 


5241 


JMP 


.-3 


5445 


1044 601 t 


TAO 


44 


5446 


7041 


CMA 


IAC 


5447 


7450 


SNA 




5452 


5256 


JMP 


G02 


5451 


3351 


OCA 


CNTR 


5452 


4337 


JMS 


GET 


5453 


4321 


JMS 


OUT 


5454 


2351 


ISZ 


CNTR 


5455 


5252 


JMP 


.-3 


5456 


1352 S02, 


TAD 


SAC 


5457 


7650 


SNA 


CLA 


5460 


5263 


JMP 


.+3 


5461 


1354 


TAO 


PERIOD 


5462 


4761 


JMS 


I 0UT1 


5463 


4337 


JMS 


GET 


5464 


4321 


JMS 


OUT 


5465 


5263 


JMP 


.-2 


5466 


1352 ERR, 


TAD 


SAC 


5467 


7700 


SMA 


CLA 


5470 


5274 


JMP 


ERGO 


5471 


1356 


TAO 


CHX 


5472 


4321 


JMS 


OUT 


5473 


5271 


JMP 


.-2 


5474 


7240 ERGO, 


CLA 


CMA 


5475 


1352 


TAD 


SAC 


5476 


3352 


OCA 


SAC 


5477 


5225 


JMP 


TRYAGN 


5500 


7200 MINS, 


CLA 




5501 


1062 


TAD 


62 


5502 


7041 


CMA 


IAC 


5503 


1352 


TAD 


SAC 


5504 


7450 


SNA 




5505 


5313 


JMP 


G03 


5506 


3351 


DCA 


CNTR 


5507 


1353 


TAO 


SPCE 


5510 


4321 


JMS 


OUT 


5511 


2351 


ISZ 


CNTR 


5512 


5307 
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PERIOD 
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OUT 
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44 
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5315 
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5545 
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GET 
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MS 
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GET 

15 
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TAD 
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-7 
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7352 

BUFFER- 1 
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5555 
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5547 
5550 
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5527 
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5537 
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